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1 . INTRODUCTION 


In  July  1973,  a FORTRAN  compiler  system^  was  announced  that 
translated  FORTRAN  IV  progrcons  into  Wang  520  calculator  code.  This 
system  was  developed  on  the  Tym-Share  Incorporated  system.  The  output 
consisted  of  a listing  of  the  translated  Wang  programs  along  with 
special  tables  needed  by  the  user  later  for  loading  the  program  into  the 
calculator. 

A new  compiler  system  is  described  here  that  contains  the  original 
system  as  a subset,  but  is  written  in  FORTRAN  IV  instead  of  SUPER 
FORTRAN.  The  following  characteristics  describe  the  changes  in  the  old 
system: 

(a)  Loader — A loading  phase  has  been  added  that  will  produce  a 
binary  card  deck  as  an  output  of  the  system.  This  deck  can  then  be 
read,  via  a mark-sense  card  reader,  into  the  Wang  520  calculator. 

(b)  Algorithm  Changes — Several  changes  that  have  been  made  in  the 
algorithms  cited  in  the  original  paper  either  remove  ambiguities  or 
improve  code  translation. 

(c)  New  Optimization — Special  attention  has  been  placed  upon  the 
availability  of  the  16  basic  registers  to  further  optimize  the  code 
space.  Also,  many  additions  have  been  made  to  the  second-pass 
optimizer,  with  special  consideration  to  the  optimization  of  array  code 
generation. 

(d)  Adaptability — Because  the  system  is  now  written  in  FORTRAN  IV, 
it  should  be  easily  adaptable  to  computers  other  than  the  IBM  360/370 
series  on  which  it  is  presently  stored. 

The  following  sections  describe  how  the  new  improved  system  is  used 
and  give  a more  detailed  account  of  the  added  characteristics. 


2.  USING  THE  SYSTEM 

The  compiler-loader  system  has  been  implemented  on  an  IBM  370/195 
computer,  accessible  at  Harry  Diamond  Laboratories  (HDL)  through  the 
IBM  1130  remote  batch  terminal.  The  following  JCL  statements  are 
necessary  to  execute  the  Wang  520  compiler-loader  system: 


Bloom  and  A.  Hausner,  FORTRAN  IV  Compiler  for  the  Wang  520 
Calculator,  Harry  Diamond  Laboratories  TM-73--15  (July  1973) 
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JOB  CARD 
ACCOUNT  CARD 

//JOBLIB  DD  DSN=C753. WANG, DISP= (OLD, KEEP) 

//  EXEC  PGM=COMPILER,REGION=280K 

//FTOlFOOl  DD  DSN=SCRATCH,DISP= (NEW, DELETE) ,UNIT=SASCR, 

//  SPACE=(1872, (133,5)) ,DCB=(RECFM=VS,LRECL=1868,BLKSIZE=1872,DSORG=DA) 
//FT06F001  DD  SYSOUT=A,DCB= (RECFM=FBA,LRECL=133 ,BLKSIZE=1596) 

//FT07F001  DD  SYSOUT=B 
//FT05F001  DD  * 

Compiler-loader  input 

/* 


The  system  output  consists  of  the  program  source  and  object  listings 
and  special  system  tables,  some  of  which  have  already  been 
described.^  If  the  loader  has  been  used,  object  programs  are  punched. 
However,  since  the  terminal  currently  in  use  receives  card  images  that 
contain  only  the  256  EBCDIC  characters,  the  punched  output  is  not  in  a 
form  directly  readable  into  the  Wang  520.  A program  has  been 
implemented  on  the  IBM  1130  to  convert  the  object  decks  to  the  Wang  520 
card  format  that  is  acceptable  to  the  mark-sense  reader.  The  following 
deck  should  be  submitted  for  local  (IBM  1130)  execution  for  obtaining 
the  card  conversion; 


Col  14  8 11  51 

//  JOB  1210  your  name 

//  BPNCH 

cards  from  loader 

blank  Wang  cards  (same  amount) 

Each  set  of  FORTRAN  and  Wang  routines  that  should  be  loaded  together 
is  called  a "load  group."  The  punched  output  of  the  compiler-loader 
consists  of  a leader  card  with  the  number  of  cards  (in  hexadecimal 
format)  in  the  object  deck  in  columns  1 and  2.  The  object  deck  follows, 
with  40  steps  punched  on  each  card,  also  in  hexadecimal  format.  The 
above  sequence  is  repeated  for  each  load  group.  Hence,  in  loading  the 
compiler-loader  system,  sufficient  blank  cards  should  be  added  behind 
the  system  deck  to  take  care  of  all  the  cards  to  be  punched. 


Bloom  and  A.  Hausner,  FORTRAN  IV  Compiler  for  the  Wang  520 
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The  Wang  cards  punched  correspond  one  to  one  to  the  cards  punched  by 
the  loader,  so  the  number  of  Wang  cards  required  is  the  same  as  the 
number  of  cards  put  into  the  program.  A maximum  of  500  cards  may  be 
converted  at  one  time. 

2.1  Compiler-Loader  Operating  System 

The  compiler-loader  represents  a very  simple  operating  system 
for  which  the  user  can  run  the  program  by  using  the  set  of  five  system 
operations  (table  I) . 

The  general  deck  setup  is 

(a)  *WANG  card(s)  for  each  external  Wang  program  included  in  a 
load  group, 

(b)  *FTC  card(s)  preceding  each  FORTRAN  program  deck  included 
in  a load  group, 

(c)  *LDR  card(s)  for  load  group, 

(d)  *END. 

(e)  Repeat  steps  (a)  to  (d)  for  any  additional  load  groups. 

(f)  *STOP  (always  the  last  card  in  the  deck) . 

FORTRAN  Program  Compilation. — Every  FORTRAN  main  program  or 
siabroutine  must  be  preceded  by  the  *FTC  control  card  (see  table  I) . The 
options  NOLIST  and  NOMAP  are  used  to  limit  the  amount  of  printed  output. 
The  option  AUTO  is  used  to  take  adveuitage  of  the  availability  of  free 
registers  in  a given  program  (see  sect.  5.1.2(d)).  The  previous  options 
need  not  be  specified. 

External  Program  Definition. — The  *WANG  card  allows  already 
existing  Wang  520  programs  to  be  combined  with  FORTRAN  programs  into  one 
final  loaded  program.  It  is  assumed  that  the  name  specified  in  the  NAME 
field  has  been  referenced  by  a FORTRAN  program.  All  the  options  are 
used  to  avoid  conflicts  with  marks  and  register  numbers  used  by  other 
programs  included  in  the  same  load  group.  Hence,  all  marks  and  nonbasic 
registers  used  by  the  external  program  must  be  included.  If  the  mark  or 
nonbasic  register  range  used  by  the  program  is  not  one  defineible 
sequence,  the  particular  option  can  be  repeated  as  often  as  necessary  to 
describe  all  the  sequences  actually  used.  If  more  than  one  card  is 
needed  for  a given  external  program,  additional  *WANG  cards  can  be  added 
with  the  NAME  field  blank. 
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TABLE  I.  OPERATING  SYSTEM  COMMANDS 


Program  name 

Operation Us^ required 

(1)  FTC  FORTRAN  program  compilation 
Options ; 

(a)  NOLIST — source  program  listing  not  printed 

(b)  NOMAP — Name  table,  translation,  and  reference 

table  not  printed 

(c)  AUTO  (nn-mm) — free  registers  in  range  nn-iran 

are  assigned  to  most  frequently  referenced 
variables  (if  nn-mm  is  not  specified  (i.e., 

AUTO) , range  01-15  is  used) 


(2)  WANG  External  program  definition  Yes 

Options  (all  must  be  specified  to  prevent 
conflicts) : 

(a)  Ennnn — mark  number  of  entry  point  to  pro- 

gram (must  be  specified) 

(b)  Snnnn — number  of  program  steps,  not 

including  END  PROG  (if  specified,  storage 
is  allocated  for  program  before  compiled 
programs) 

(c)  Mnnnn-mmmm — Mnnnn  is  mark  number  used  by  pro- 

gram (if  -mmmm  is  specified,  entire  range  is 
excluded  by  loader) 

(d)  Rnnn-mmm — Rnnn  is  nonbasic  register  number 

used  by  program  (if-mmm  is  specified, 
entire  range  is  excluded  by  loader) 


(3)  LDR  Assigns  storage  and  marks  to  compiled  programs  and  Optional 

satisfies  external  references  between  programs 

Options: 

(a)  Mnnnn — starting  mark  number  (if  not  specified, 

first  available  mark  is  used) 

(b)  Rnnn — starting  nonbasic  register  number  (if  not 

specified,  first  available  number  is  used) 

(c)  Ennn — entry  point  mark  number  (if  not  specified, 

first  available  entry  point  number  is  used) 


(4) 

END 

End 

of 

load  group 

No 

(5) 

STOP 

End 

of 

compiler  input 

No 

General  Form 

Co)  1 Col  2-5  Col  8-13  Col  16-71 

* OP  NAME  OPT, , 0PT_ , . . . , OPT 

12  n 

where  OP  is  the  operation,  NAME  is  a prograun  name,  and  OPT^  are  option  fields,  which 
may  appear  in  any  order. 
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Program  Loading. — The  *LDR  card  is  used  to  assign  storage  and 
marks  to  compiler  programs  and  satisfy  external  references  between 
programs.  If  NAME  is  specified,  it  must  have  appeared  previously  on  an 
*FTC  card.  If  NAME  is  blank,  all  programs  compiled  or  declared  external 
since  the  last  *END  Ccird  are  loaded.  If  one  of  a load  group  is  neuned  on 
an  *LDR  card,  each  program  must  have  a separate  *LDR  card  with  the  name 
used  on  the  corresponding  *FTC  card.  None  of  the  options  need  be 
specified. 


End  of  a Load  Group. — The  *END  card  marks  the  end  of  a group  of 
progreUtis  to  be  loaded  together.  Whenever  the  loader  is  used,  an  *END 
Ccurd  must  follow  the  *LDR  card(s),  to  complete  the  loading  process. 
There  may  be  several  load  groups  in  one  job,  each  terminated  by  an  *END 
card. 


End  of  Compiler  Input. — The  *STOP  card  is  used  at  the  end  of 
all  compiler-loader  input  to  halt  execution  of  the  system. 

2.2  Control  Card  Examples 

(a)  Compile  a single  program  and  do  not  load.  All  listings  are 

desired. 

*FTC  MAIN 

(FORTRAN  deck) 

★END 

★STOP 

(b)  Compile  emd  load  subroutines  A and  B.  Do  not  get  listings 
of  the  trcuislation;  specify  automatic  register  allocation  for  A.  In  the 
same  job,  compile  emd  load  subroutine  C with  external  progreun  D. 
Program  D uses  marks  0003  and  0010-0103  and  registers  016-047 ; its  entry 
point  is  1010,  auid  it  uses  156  steps. 

*FTC  A NOMAP,  AUTO 

(subroutine  A deck) 

★FTC  B NOMAP 

(subroutine  B deck) 

★LDR 

★end 

*WANG  D ElOlO,  S0156,  M0003,  MOOlO-0103,  R016-047 

★FTC  C 

(subroutine  C deck) 

*LDR 

★END 

★STOP 
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(c)  Compile  and  load  program  INTG.  Registers  03-07  are  to  be 
used  for  automatic  register  allocation.  No  source  listing  is  desired. 
The  loaded  program  uses  marks  beginning  with  1200  and  registers 
beginning  with  032.  The  entry  point  for  INTG  is  1104. 

*FTC  INTG  AUTO (03-07),  NOLIST 

(FORTRAN  deck) 

*LDR  M1200,  R032,  E1104 

*END 

*ST0P 

(d)  Compile  and  load  programs  FTC  and  DERIV  together. 
Program  FTC  should  use  marks  beginning  at  0800  and  registers  beginning 
at  025.  Program  DERIV  uses  marks  starting  at  0900  and  registers 
beginning  at  025.  Even  though  the  starting  register  number  is  the  same, 
the  loader  does  not  assign  the  same  registers  to  both  programs.  The 
first  register  used  by  DERIV  is  the  first  register  after  25  not  used  by 
FCT. 

*FTC  FCT 

(FORTRAN  deck) 

*FTC  DERIV 

(FORTRAN  deck) 

*LDR  FCT  M0800,  R025 

*LDR  DERIV  M0900,  R025 

*END 
*STOP 

2 . 3 FORTRAN  Implementation  and  Limitation 

The  FORTRAN  IV  siibset  implemented  by  the  Wang  compiler-loader 
system  has  been  completely  described.^  To  make  this  report  as 
self-sufficient  as  possible,  a general  list  of  unimplemented  FORTRAN  IV 
capabilities  is  reported  here  (table  II) . 


Bloom  and  A,  Hausner , FORTRAN  TV  Compiler  for  the  Wang  520 
Calculator f Harry  Diamond  Laboratories  TM-73-15  (July  1973) , section  2. 
Changes  in  the  restrictions  are  noted  in  section  4 of  TM-73-15S. 
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TABLE  II.  FORTRAN  PROGRAM  STATEMENTS  REQUIRINT  EDITING^ 


Subroutines 
not  provided 

Not  implemented 
and  ignored 

Not  implemented 
with  error  result 

Must  be  followed  by 
at  least  one  blank 

OVERFL ( J) 

ENDFILE 

ASSIGN  N TO  I 

DIMENSION 

DVCHK(J) 

PUNCH 

BLOCK  DATA 

COMMON 

SSWTCHd.J) 

PRINT 

COMPLEX 

INTEGER 

SLITET(I) 

PJEWIND 

DATA 

REAL 

ERF(X) 

NAMELIST 

ENTRY 

DOUBLE 

GAMMA (X) 

BACKSPACE 

CALL  EXIT 

PRECISION 

ALGA>yv(X) 

FORMAT 

EXTERNAL 

LOGICAL 

CEXP(X) 

GOTO  I,  (Nl, . . . , 

SUBROUTINE 

CLCX3(X) 

Machine-dependent 

NM) 

CALL 

CSIN(X) 

functions  ignored 

RETURN  1(1  not 

FUNCTION 

CCOS(X) 

AND(X,Y) 

blank) 

CABS (X) 

OR(X,Y) 

Arithmetic  State- 

CSORT(X) 
CMPLX(X1,X2) 
CONJG (X) 
AIMAG(X) 
REAL(X) 
DUMP(A  ,B 
PDUMP  ^ ^ 

(A^ , , . • . ) 

COMPL(X) 
BOOL (X) 

ment  Function 

Special  cases; 

(a)  Function  names  not  allowed  as  arguments  of  a subroutine  or  function 
subprogram.  Subscripted  variables  allowed  when  not  outputs.  Minimum  of 
1 and  maximum  of  15  arguments.  Function  nesting  limited  to  5. 

(b)  Values  of  DO  indices  are  not  available  outside  of  the  loop. 

(c)  Varicibles,  constants,  reserve  words,  or  special  operators  cannot  be 
continued  on  the  next  line. 


^From  H,  Bloom  and  A,  Hausner , Harry  Diamond  Laboratories  TM-73-15, 


In  general^  all  the  important  facilities  in  FORTRAN  IV  have 
been  implemented;  perhaps  the  only  two  deficiencies  are  no  complex 
arithmetic  and  only  single-dimension  arrays.  However^  these 
restrictions  are  sensible  when  the  computer  being  used  has  the  limited 
size  of  a Wang  520.  Table  III  summarizes  the  specific  restrictions. 
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TABLE  III.  SET  OF  FORTRAN  IV  CAPABILITY  RESTRICTIONS* 


(1)  Complex  constants  are  not  allowed. 

(2)  Arrays  can  be  only  one  dimensional. 

(3)  Trigonometric  functions  arguments  such  as  SIN  have  a magnitude 
restriction  of  10  radians. 

(4)  Computed  GOTO  has  a limit  of  20  statement  numbers.  If  the 
index  is  out  of  range,  the  default  is  to  the  last  number. 

(5)  The  DO  index  value  is  not  available  outside  the  DC  loop.  The 
maximum  DO  nest  level  is  5 . 

(6)  The  "n"  in  PAUSE  or  STOP  may  be  only  a single  decimal  digit. 

(7)  Maximum  and  minimum  built-in  functions  can  have  two  arguments 
only. 

(8)  User-defined  subprograms  do  not  supersede  system  routine 
names . 

(9)  Variable  dimensions  are  not  allowed. 

(10)  Only  one  labeled  COMMON  can  appear  on  a given  C0M310N  card. 


^Summarized  from  H.  Bloom  and  A,  Jiausner , Harry  Diamond  Laboratories 
TM-7 3-1 5 , section  2. 

The  compiler-loader  system  assumes  that  the  FORTRAN  decks 
submitted  are  free  of  syntax  errors  and  thus  does  not  perform  any  degree 
of  syntax  checking  itself.  Therefore,  it  is  possible  for  a run  to  fail, 
and  the  computer  will  list  an  error  that  is  completely  meaningless  to 
the  user.  In  this  case,  the  user  should  satisfy  himself  that  his  decks 
are  running  programs. 

2.4  User-Supplied  Compiler  Commands 

As  described  in  detail  previously,^  the  system  allows  the  user 
to  include  special  compiler  commands  (table  IV)  to  improve  his  code,  if 
he  desires.  These  commands  appear  directly  in  the  FORTRAN  deck  and  are 
entered  as  special  comments.  The  commands  T,  R,  and  S are  used  for 
register  optimization.  Since  the  automatic  register  assigner  sxibsystem 
(AUTO’,  sect.  3.3)  was  added,  there  is  little  need  for  the  user  to  assign 
variables  to  basic  registers.  However,  the  option  is  available  and  does 


Bloom  and  A.  Hausner,  FORTRAN  IV  Compiler  for  the  Wang  520 
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TABLE  IV.  USER-SUPPLIED  COMPILER  COMMANDS* 


1. 


Command 


Effect 


T Z 

(Top  assigning) 


The  top  register  in  the 
from  15  to  Z (<15)  . 


program 


is 


changed 


2.  R N(V^, . .. ,V  ) 
(Equivalencing) 


The  FORTRAN  variables  all 

assigned  the  same  basic  register  N(00<.Nj<15)  . 


3.  S name(C  , . . . ,C  ) , 
d,  R1-R2,  $Z 
(Subprogram 
specification) 


(a)  "name"  is  the  name  of  a subprogram  to  be 

called  as  name (a^, . . . ,a  ). 

0 m 

(b)  is  a special  symbol  that  determines 
that  property  of  the  corresponding  actual 
subprogram  call  argument  as  follows : 

C^=B; argument  a.  is  both  an  input  and 
output  variable  (a.  must  be 
nonsubscripted) . ^ 


I ; argument 
only. 


a. 


1 


is  input  expression 


O;  (letter  0)  argument  a.  is  output 
variable  only.  (a^  must  be 

nonsubscripted) . ^ 


E; argument  a.  is  empty  variable,  i.e., 
neither  input  nor  output. 

(c)  d is  the  maximum  number  of  nested  DC 
indices  in  "name."  If  d is  omitted,  all 
necessary  DO  indices  of  the  calling 
program  will  be  saved  and  restored. 

(d)  R1-R2  is  the  range  of  basic  registers  to 
be  saved  in  the  calling  program.  If 
Rl=R2=0,  no  registers  are  saved.  If  the 
range  is  omitted,  all  variables  assigned 
to  basic  registers  in  the  calling  program 
will  be  saved  and  restored. 


(e)  Z is  the  top  register  specified  in  a T 
statement  in  "name."  If  $Z  is  omitted, 
Z=15. 


(f)  d,  R1-R2,  and  $Z  fields  may  be  listed  in 
any  order . 


^From  H.  Bloom  and  A.  Hausner,  Harry  Diamond  Laboratories  TM-^73-15 
(July  1973) 
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TABLE  IV.  USER-SUPPLIED  COMPILER  COMMANDS^  (CONT'D) 


4. 


Command 


W V (format,  N,  M)  (a) 

(I/O  specification) 

(b) 


Effect 


V is  the  name  of  a FORTRAN  variable. 

Format  is  any  four  digit  allowable 
format  that  can  follow  the  PRINT 
(table  33  of  HDL-TM-73-15) . 


Wang 

step 


(c)  N is  the  print-on-read  indicator. 


N = 0;  do  not  print  input  on  P^AD. 

N = 1;  print  input  on  READ. 

(d)  M is  the  spacing  indicator. 

M = 0;  no  spaces. 

M - 1;  space  before  printout. 

M = 2;  space  after  printout. 

M = 3;  space  before  and  after  printout. 


(e)  If  N = 0,  M is  ignored  on  a READ. 


(f)  If  format  is  0015,  N and  M are  used  with 
the  previously  defined  default  format 
(table  19  of  HDL-TM-73-15). 


(a)  All  commands  must  start  with  COMPILE  in  columns  1-7. 

(b)  Card  commands  cannot  be  continued  on  next  line. 

Cc)  Only  one  S command  can  appear  for  each  "name." 

(d)  Any  FORTRAN  variable  can  appear  in  at  most  one  R command. 

(e)  Subprogram  arguments  cannot  appear  in  an  R command. 

(f)  Only  one  register  can  be  specified  in  one  R command. 


^From  H.  Bloom  and  A,  Hausner,  Harry  Diamond  Laboratories  TM-73-15 
(July  1973) 
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allow  for  the  most  efficiently  generated  code.  The  command  T is  usually 
used  with  respect  to  the  command  S,  in  order  to  line  up  arguments  in  the 
called  routine  with  the  argument  register  storage  in  the  calling 
routine.  The  S command  is  mainly  designed  to  minimize  the  amount  of 
code  that  must  be  generated  each  time  a call  is  translated.  Examples  of 
ways  that  the  R,  S,  and  T commands  can  be  used  have  been  given  earlier. 

Since  format  statements  are  not  implemented  in  this  system,  the 
W command  allows  the  user  to  make  full  use  of  the  I/O  capabilities  of 
the  Wang  520  printer.  Hence,  whenever  a variable  is  referenced  in  an 
I/O  statement,  the  indicated  W format  is  used  in  the  translation,  if  it 
exists. 


3.  DESCRIPTION  OF  SYSTEM  STRUCTURE 

The  compiler  is  divided  into  six  subsystems:  (1)  encoder,  (2) 

parser,  (3)  automatic  register  assigner,  (4)  translator,  (5)  optimizer, 
and  (6)  loader.  A brief  description  of  each  system  is  given  below. 

3 . 1 Encoder 


The  encoder  takes  the  FORTRAN  program  string  code  and  generates 
an  encoded  program  string.  It  performs  the  following  special  features: 


(a)  Generates  special  operator  code  (for  example,  X^,  1/X,  II, 


power  of  10  shifts)  to  take  advantage  of  the  Wang  operator  set. 


i 


(b)  Sets  up  COMPILE  option  tables. 


(c)  Eliminates  FORTRAN  formats  from  code  and  removes  format 
number  and  device  index  from  I/O  statement. 


(d)  Removes  commands  not  allowed  in  version. 

(e)  Generates  special  end-of-DO  (CONTINUE)  statement  with  new 
label  that  also  replaces  the  DO  statement  label.  This  is  a new  feature 
that  eliminates  ambiguities  in  the  pretranslation  of  a FORTRAN  statement 
into  a multiline  code  (i.e.,  for  I/O,  certain  LOGICAL  IF  and  function 
calls) . 

3.2  Parser 


The  parser  talces  the  encoded  string  and  generates  a reverse 
Polish  string.  It  performs  the  following  special  features; 


%.  Bloom  and  A.  Hausner,  FORTRAN  IV  Compiler  for  the  Wang  520 
Calculator,  Harry  Diamond  Laboratories  TM-73-15  (July  1973),  table  12, 
section  3.2.1. 


15 


(a)  Reorders  expressions  having  a commutative  binary  operator 
(+^  .AND.^  .EQ.^  .OR.,  .NE.),  if  there  is  a reduction  in  the  number 
of  registers  needed  to  compute  the  expression.  The  new  algorithm 
replaces  the  one  mentioned  previously^  in  that  it  simply  checks  the 
expression  with  the  arguments  in  either  position  to  determine  the  best 
ordering . 


(b)  Reorders  commutative  expressions,  provided  there  is  no 
change  in  the  number  of  registers  under  the  following  conditions  (in 
order  of  priority) : 


(1)  Reorders  the  assignment  expression  so  that  a variable 
assigned  in  the  previous  statement  is  closest  to  the  left  of  the  present 
expression.  Doing  so  allows  the  recall  suppression  feature  to  be  more 
effective;  for  example,  statements  k = l,  I = J + K will  be 
pretranslated  asK=l,  I=K+J. 


(2)  Reorders 
consecutive  operators 
B becomes  X = B - A. 


unary 


the  expression 

+ with  . 

binary 


to  replace  the  two 
for  example,  X = - A + 


(c)  Defines  code  for  (F,  (A,  and  (M  to  indicate  the  end  of  a 

function,  array,  and  macro  argument  list,  respectively. 

(d)  Checks  LOGICAL  IF  statement  for  function  or  I/O  statement 
expansion  candidates.  If  expansion  is  required,  the  statement  is 
altered.  For  example, 

IF(X.  EQ.l)  Y(I)  = F(X) 

becomes 

IF(X.  EQ.l)  GOTO  LI 
GO  TO  L2 
LI  Y(I)  = F(X) 

L2  next  statement 

where 


LI  and  L2  are  system-generated  labels. 


Bloom  and  A.  HausneTf  FORTRAN  IV  Compiler  for  the  Wang  520 
Calculator f Harry  Diamond  Laboratories  TM-73-15  (July  1973),  table  12, 
section  3.2.1. 


16 


(e)  Stores  equivalence  and  common  statements  into  system 

tables. 

(f)  Checks  for  function  or  I/O  expansion.  If  the  statement  has 
a labels  the  label  goes  with  the  first  statement  in  the  expansion.  For 
example , 

10  Y(I)  = F(X)  becomes  10  T.OO  = F(X) 

Y(I)  = T.OO 

10  READ  ( ) X,  Y becomes  10  READ  { ) X 

READ  ( ) Y 

(g)  Defines  the  array  repeat  operator  A^^  if  assignment  array 
appears  on  the  right  side  of  the  assignment.  For  example, 

X(I)  = X(I)  + B/X(I) 


becomes 


X(I)  = \ ^ B/A^ 

(h)  Looks  for  the  pattern  "expression  * expression"  (appearing 
in  FORTRAN  format)  and  replaces  it  with  (expression) **2. 

(i)  Looks  for  the  pattern  ++  or  **  (appearing  in  reverse 
Polish).  If  it  is  found,  the  second  operand  of  the  first  operator  is 
switched  with  the  operator.  For  example, 

the  FORTRAN  (A  + B)  + (C*  D) , first  in  Polish  (AB+CD*+) , 
and  then  (CD*AB++)  after  reordering  becomes 


the  FORTRAN  ( (C  * D)  + A)  + B,  or  in  Polish  (CD*A+B+) . 

3.3  Automatic  Register  Assignor  (AUTO) 

The  AUTO  system  assigns  auxiliary  registers  (i.e.,  work 
registers,  DO  index  registers,  and  register  00)  to  common  subexpressions 
and  temporary  variables  and  free  registers  to  other  variables  if  so 
specified.  The  discussion  of  register  optimization  is  very  important 
and  is  handled  specially  in  section  5.1.  This  system  is  new  to  the 
translator,  but  the  system  can  greatly  reduce  the  program  size  by  trying 
to  optimize  the  use  of  the  basic  registers. 


17 


3.4  Translator 


parsed  string  and  generates  the 
many  local  optimization  features 


The  translator  takes  the 
relocatable  Wang  code.  There  are 
performed  at  this  point. ^ 

The  techniques  used  in  translating  each  type  of  FORTRAN 
statement  into  corresponding  Wang  Code  have  been  described^  in 
detail.  Even  if  the  compiler-loader  system  is  not  used,  it  is 
worthwhile  to  learn  the  translation  algorithms  and  apply  them  when 
writing  the  Wang  code,  especially  in  the  case  of  conditional  transfers 
and  DO  loops.  The  original  report  also  describes  the  basic  translator 
scheme  for  setting  up  the  calculator  as  a computer  structure  and  use  of 
the  registers  in  implementing  the  translation  algorithms. 


3.5  Optimizer 

The  optimizer  (or  second-pass  optimizer)  takes  the  relocatable 
Wang  code  generated  by  the  translator  and  searches  for  certain  patterns, 
in  order  to  reduce  the  number  of  steps.  The  seven  passes  made  through 
the  code  have  also  been  discussed.^  The  current  version  of  the 
optimizer  makes  an  additional  eight  passes  to  try  to  further  reduce  the 
number  of  steps.  These  are  fully  discussed  in  section  5.2. 


3.6  Loader 


The  loader  combines  the  whole  relocatable  Wang  code  generated 
from  the  FORTRAN  programs  with  any  externally  defined  Wang  programs  and 
stores  the  whole  code  together  in  absolute  form  by  use  of  the  following 
order : 


(a)  External  Wang  programs 

(b)  FORTRAN  programs 

(c)  Large  number  storage  in  nonbasic  registers 

(d)  Nonbasic  variable  register  storage. 

The  loader  outputs  a set  of  punched  cards  representing  the 
object  deck  of  a loaded  program  module  (sect.  2) . 


Bloom  and  A.  Hausner,  FORTRAN  IV  Compiler  for  the  Wang  520 
Calculator,  Harry  Diamond  Laboratories  TM-73-15  (July  1973),  tables  15, 
17,  18. 
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An  optimization  algorithm  within  the  loader  attempts  to  use  the 
same  nonbasic  registers  for  variables  in  different  programs,  if  doing  so 
is  possible.  Variables  (including  arrays)  that  are  candidates  for 
sharing  a nonbasic  register  with  a variable  (or  array)  in  another 
program  must  obey  the  following  property.  Between  the  first  and  last 
reference  to  the  variable,  there  can  be  no  call  to  another  subprogram. 
This  restriction  is  necessary  to  keep  the  value  of  the  assigned  nonbasic 
register  from  changing  in  the  called  subprogram  and,  hence,  storing  an 
incorrect  value  for  the  variable  used  in  the  calling  program. 


4.  CHANGES  IN  ALGORITHMS 


Aside  from  many  new  features  in  the  compiler  system,  there  have  been 
many  changes  also  in  the  algorithms  already  developed  and  described.^ 

(a)  The  code  generated  from  the  nonbasic  operators  .LE.,  .LT., 
.GE.,  and  .GT.^  has  been  shortened,  and  the  code  for  the  .EQ.  and  .NE. 
simple  LOGICAL  IF  statement  has  been  changed,  to  avoid  the  use  of 
register  00.  Table  V shows  the  changes. 

TABLE  V.  CHANGES  IN  LOGICAL  CODE 


.LE. 

.LT. 

.GE. 

.GT. 

IF  (J.EQ.K)  S 

IF  (J.NE.K)  S 

- L 

- L 

“ L 

- L 

recall  J 

recall  J 

T L 

T L 

E 1 

E 1 . 

ST  L 

ST  L 

J IF  0 

J IF  0 

SP-ST  L 

SP-ST  L 

recall  K 

recall  K 

E 0 

J IF  + 

J IF  0 

J IF  0 

- L 

- L 

J IF  + 

E 1 

J IF  + 

J IF  + 

J NE  0 

J IF  0 

E I 

ST  L 

T L 

GO 

translate  S 

translate  S 

ST  L 

RE  L 

GO 

T L 

RE  L 

RE  L 

RE  L 

(b)  The  code  used  for  generating  function  expansion  ^ has  been 
altered  so  that  if  a label  appears  in  the  original  statement,  this  label 
eventually  appears  in  the  first  statement  of  the  expansion. 
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(c)  The  code  for  generating  DO  loops ^ has  been  altered.  The 
special  DO  label  Lnn  is  now  used  as  the  label  on  the  last  statement  of 
the  DO  loop.  (A  special  CONTINUE  statement  is  added,  see  sect.  3.1(e).) 
The  combination  of  changes  (b)  and  (c)  removed  the  problem  of  not 
allowing  function  calls  to  be  in  the  last  statement  of  a DO  loop. 

(d)  The  special  DO  label  name  generated  is  now  D.nm. 

(e)  The  system  label  generated  is  now  LABnm. 

(f)  The  system-created  variable  (ZZZZnm)  is  now  T.nm. 

(g)  The  subprogram  save  variable  (AAAAnm)  is  now  AAA.nm. 

5.  NEW  OPTIMIZATION  FEATURES 

Since  the  development  of  the  original  compiler,  the  major  change  in 
the  new  system  (besides  the  loader  phase)  has  been  the  concern  for  more 
step  optimization.  The  work  lies  in  two  major  areas: 

(a)  maximvim  use  of  the  basic  registers 

(b)  recognition  of  repeatable  patterns  in  the  code  so  that  copies 
of  the  patterns  can  be  eliminated  (e.g.,  generation  of  the  same  array 
element  in  two  statements) . The  two  types  of  optimization  take  place  in 
the  AUTO  and  optimizer  subsystems,  respectively. 

The  user  previously  had  to  decide  whether  or  not  he  needed 
optimization  tricks.  They  are  now  automatically  performed  by  the 
compiler  (table  VI) . 


%.  Bloom  and  A.  Hausner,  FORTRAN  IV  Compiler  for  the  Wang  520 
Calculator,  Harry  Diamond  Laboratories  TM-73-15  (July  1973). 
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TABLE  VI.  ADOPTION  OF  OPTIMIZATION  TRICKS 


Type 

Section  with 
reference  ^ 

Using  the  basic  registers 
effectively 

5.1  (6.2.1) 

Arithemetic  IF  statements 

5.2  (6.2.2) 

DO  index  suppression 

3.2-b-l  (6. 2.4-2) 

Low  index  (of  DO)  is  an 
expression 

5.1.2-b  Case  2 (6. 2. 4-3) 

Register  00  is  not  used 
(upper  index  is  expression) 

5.1.2-b  Case  3 (6. 2. 4-4) 

Recall  array  suppression 

5.1.2-a  (6. 2. 5-1) 

Transposing  terms  for  recall 
suppression 

3.2-b-l  (6. 2. 6-1) 

Transposing  terms  to  lower 
work  register  requirements 

3.2-a,  3.2-i  (6. 2.6^3) 

Trailing  zeros 

5.2  (6.2. 7-1) 

Constant  set 

5.2  (6.2. 7-2) 

Decimal  suppression 

5.2  (6.2. 7-4) 

Squaring 

3.2-h  (6. 2. 8-3) 

ZZZZnn  equivalencing 

5.1.2-c  (6. 2.8-4) 

^Section  number  in  parentheses  is  from  HDL-TM-73-15 . 

5.1  Register  Optimization 

As  mentioned  in  section  3.3,  the  most  important  resources  in 
the  calculator  are  the  16  basic  registers  that  play  an  important  role  in 
the  statement  translation.  Some  of  the  optimization  concepts  were 
obtained  from  Cries, ^ Hopgood,^  and  Rustin.**  However,  it  is  important 
also  that  some  of  the  registers  be  used  for  storing  the  values  of 
variables.  There  are  three  very  critical  savings: 

^D.  GrieSf  Compiler  Construction  for  Digital  Computers,  New  York: 
John  Wiley  and  Sons  (1971) . 

^F.  R.  A.  Hopgood,  Compiling  Techniques,  New  York:  American  Elsevier 
Publishing  Co.  (1969). 

^R.  Rustin,  ed..  Design  and  Optimization  of  Compilers,  Englewood 
Cliffs,  NJ:  Prentice  Hall,  Inc.  (1972). 
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(a)  Since  "nonbasic"  registers  are  constructed  by  the 
concatenation  of  eight  program  steps,  these  steps  can  be  saved  if  the 
variables  can  be  stored  in  basic  registers,  instead. 

(b)  Many  of  the  step  optimization  features  (such  as  simple 
updating)  can  be  applied  only  if  the  varible  is  in  a basic  register. 

(c)  Execution  time  is  saved,  not  only  from  the  above  types  of 
step  savings,  but  also  by  the  use  of  basic  register  store-recall  rather 
than  "nonbasic**  store-recall. 

The  goal  of  register  optimization  is  to  make  use  of  the  basic 
registers  when  they  are  not  being  used  in  the  ordinary  statement 
translation.  Let 

M = the  maximum  number  of  work  registers  used  in  any  one 
statement , 

M-  = the  maximum  number  of  DO  index  registers  needed  at  any 
^ . d 

time, 

LT  = the  first  available  work  register, 

= the  number  of  work  registers  needed  for  statement  i, 

D.  = the  DO  index  level  at  statement  i. 

1 

For  statement  i,  the  registers  LT-W.  through  LT-M  +1  and 
through  D.+l  are  free  for  some  other  use.^  In  addition,  register  00 
might  also  be  free.  Hence,  these  registers  can  be  used  for  storing  the 
values  of  variables  at  selected  points  in  the  program.  In  programs 
where  there  is  a large  value  for  W.  (such  as  4)  or  Dj^  (such  as  3)  , 
register  optimization  can  result  in  tremendous  savings  of  space,  since 
these  peak  values  occur  infrequently  in  the  program  (for  W. , usually 
once  or  twice) . However,  care  must  be  taken  in  the  choice  of  which 
varicJDles  to  use  with  which  registers.  The  next  section  lists 
definitions  of  terms  needed  in  the  description  of  the  register 

optimization  that  is  performed  by  the  compiler. 

5.1.1  Definitions 

(a)  Auxilliary  register — any  register  that  at  any  time  was 
used  as  either  a work  register  (i.e.,  for  computing  expressions)  or  a DO 
index  register.  Register  00  is  also  considered  auxiliary.  These 
registers  are  used  at  least  once  in  the  ordinary  statement  translation 
described  in  section  3.4.  An  auxiliary  register  is  considered  available 
in  a given  statement  span  if  it  has  not  been  previously  used  in  any 
statement  within  the  span. 
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(b)  Free  register — free  registers  F are  those  defined  by 
LT-M^  > > M^.  These  registers  are  never  used  in  the  ordinary 
statement  translation;  hence,  they  are  entirely  free  to  be  used  for 
variable  register  assignment. 

(c)  Variable  domain — the  span  of  statements  from  the  first  to 
the  last  referencing  the  variable  in  the  program.  If  the  variable 
appears  outside  a DO  loop  and  its  last  reference  is  within  the  DO  loop, 
its  domain  is  extended  to  include  the  entire  DO  loop. 


(d)  Common  subexpression  domain — the  span  of  statements  from 
the  first  to  the  last  that  includes  the  subexpression.  Within  the  span, 
the  following  criteria  must  be  satisfied: 


case: 


(1)  No  statement  labels  except  for  the  following  special 


IF  (arithmetic  expression)  LI,  L2,  L3 
statement  where  i = 1,  2,  or  3. 

(This  criterion  is  somewhat  restrictive,  but  greatly 
simplifies  the  pattern  searching.  Also,  it  eliminates  the  need  to  move 
the  expression  backwards  in  the  program  to  make  certain  that  no  j\imp  can 
occur  over  the  first  appearance  of  the  expression.) 


expression 


(2)  No  new  assignment  involving  a variable  used  in  the 


(3)  No  subprogram  references. 

(e)  Closed  block — a span  of  statements  obeying  the  following 

properties: 

(1)  No  backward  transfers 

(2)  No  transfers  from  outside  the  block  to  within  the 

block 

(3)  DO  loops  may  appear  within  a closed  block 

(4)  No  function  or  svibroutine  call  can  appear  except  as 
the  first  statement  of  the  closed  block. 
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(without  properties  (1)  and  (2),  the  span  of  the  closed 
block  would  be  inefficiently  expanded  to  include  the  whole  range  from 
the  transfer  back  to  the  label  reference.  Doing  so  would  eliminate  most 
of  the  register  optimization.  Property  (3)  is  included  because  of  the 
popular  usage  of  the  DO  loop,  which  usually  dominates  every  segment  of 
any  progr2on.  Property  (4)  is  necessary  because  the  compiler  does  not 
know  what  registers  may  be  used  in  the  subprogram;  hence,  it  would  have 
to  save  all  registers  used  up  to  that  point  and  result  in  extra  steps 
generated  by  the  save  code . ) 

(f)  Local  domain  of  a variable — any  span  of  statements  that 
references  the  variable  and  obeys  the  following  properties: 

(1)  The  first  statement  is  an  assignment  to  the  variable 
(not  as  a result  of  a LOGICAL  IF) . This  cannot  be  an  update  (Y  = f (Y) ) . 

(2)  All  properties  of  a closed  block  must  be  obeyed. 

(3)  The  span  includes  all  references  after  the  first 
assignment  that  obey  properties  (1)  and  (2) . 

(The  properties  define  a region  in  which  a variable  is 
totally  defined,  since  the  first  statement  must  be  an  absolute 
assignment. ) 

(g)  Temporary  varicible — a variable  whose  domain  lies  entirely 
within  a single  closed  block.  In  addition,  the  variable  may  not  have 
appeared  in  a COMMON,  FUNCTION,  SUBROUTINE,  or  EQUIVALENCE  statement, 
nor  be  an  array  name.  A good  example  of  a temporary  .variable  is  the 
switch  variable  used  for  arrays: 

TEMP  = A(I),  A(I)  = A(I+1),  A(I+1)  = TEMP. 

The  variable  TEMP  has  no  real  importance,  except  to  be  used  as  temporary 
storage.  Certainly,  TEMP  should  be  stored  in  a basic  register,  if 
possible. 


(h)  Global  variable — a variable  whose  domain  includes  one  or 
more  closed  blocks.  These  variables  are  usually  heavily  referenced 
throughout  the  program;  hence,  they  cannot  fit  within  any  defined 
domain.  However,  it  would  be  especially  beneficial  if  basic  registers 
could  be  assigned  them  because  of  the  step  savings  in  the  referencing. 

(i)  Local  variable — A variable  over  whose  domain  exists  one 
or  more  local  domains.  The  local  variable  obeys  the  properties  of  a 
temporary  variable.  Local  variables  are,  in  most  cases,  special 
examples  of  temporary  variables.  In  general,  the  programmer  has  used 
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the  name  of  a variable  that  is  constantly  used  for  temporary  storage  in 
many  places  in  the  program.  The  local  domain  definition  implies  that 
the  variable  is  entirely  defined  within  a small  local  region^  just  as  in 
the  example  given  for  the  temporary  variable — see  (g)  above.  This  local 
definition  may  occur  several  times  for  the  same  variable  in  one  progrcim. 

5.1.2  Register  Assignment 

As  stated  in  section  3.3^  auxiliary  registers  are  assigned  by 
the  system  in  its  process  of  translating  the  FORTRAN  statements.  This 
section  is  concerned  with  the  assignment  of  registers  for  varicQ^les. 
The  system  does  assign  registers  to  all  arguments  of  a subprogram 
(fig.  1).  It  can  allow  the  programmer  to  specifically  assign  variables 
to  registers  under  program  control  (sect.  2.4).  The  system  register 
assignments  are  made  in  the  following  order: 


00 

01 


(L) 


LT- 


Register  (basic)  number 


Register  (nonbasic)  number 


T— 15 


SPECIAL  WORK  REGISTER 


DO  INDEX  REGISTERS 


WORK  REGISTERS 


ARGUMENTS 
OF  SUBROUTINES 


16- 


0 

step  No. 


DATA 

MEMORY 


PROGRAM 

STORAGE 


Figure  1.  Configuration  of  translator  scheme. 
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(a)  Assignment  of  common  siibexpressions — If  auxiliary 
registers  are  available,  common  subexpressions  are  assigned  to  them  in 
the  following  order.  Let  p = the  number  of  occurrences  times  the 
number  of  tokens  (e.g.,  operators  and  operands)  in  an  expression. 
Hence,  p is  a rough  estimate  of  the  number  of  steps  that  can  be  saved  by 
replacing  all  occurrences  of  a common  subexpression  with  a variable 
reference  and  initially  assigning  the  variable  the  value  of  the 
expression.  The  expressions  are  considered  in  decreasing  values  of  p. 
The  register  used  is  no  longer  available  in  the  entire  expression 
domain.  If  a common  subexpression  cannot  be  assigned  a register,  it  is 
ignored . 


Common  subexpression  register  assignment  was  chosen  first 
because  the  domain  is  so  restrictive  and  huge  savings  of  steps  are 
possible.  Experience  has  shown  that  most  subexpressions  are  array 
elements  (such  as  A(I+1)).  In  this  compiler,  it  takes  at  least  six 
steps  to  generate  an  array  element.  Hence,  if  the  element  appears  three 
times,  at  least  9 steps  can  be  saved  (12  steps,  less  two  extra  recalls 
and  one  store) . Expressions  are  ignored  if  there  are  no  free  auxiliary 
registers,  since  use  of  a temporary  variable  for  storage  requires  the 
generation  of  a nonbasic  register  that  requires  an  additional  eight 
steps.  Experience  has  shown,  however,  that  a free  auxiliary  register  is 
almost  always  available. 

(b)  Assignment  of  a temporary  variable — temporary  variables 
are  assigned  to  auxiliary  registers,  if  available,  in  the  order  of 
decreasing  frequency  of  references.  The  register  used  is  no  longer 
available  in  the  entire  variable  domain. 

Approximately  20  percent  of  the  variables  in  an  average 
program  obey  the  definition  of  a temporary  variable.  However,  only  a 
few  of  these  variables  have  an  auxiliary  register  available  over  their 
domain,  unless  the  domain  is  relatively  small  in  its  span. 

In  three  special  cases,  the  auxiliary  register 
availability  definition  can  be  overridden  and  provide  even  more 
optimization.  The  cases  follow: 

Properties  of  Variables  Considered  Special  Cases 

(1)  Exactly  two  references  are  in  consecutive  statements 
(in  the  entire  program) . 

(2)  First  reference  is  an  assignment. 

(3)  Second  reference  is  not  in  a labelled  statement. 
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Cases 


or 


(1)  Strictly  temporary  variable: 


X=.  . 

A=X. 


Register  LT  is  used  for  X if  A is  nonarray 
and  LT-1  is  used  if  A is  an  array. 


K=.  . . I 

> Register  LT  is  used  for  K. 

A(K)=.  . .) 

For  this  case,  X (or  K)  would  have  been  recalled  and 
then  stored  in  LT  anyway;  instead,  the  values  can  be  stored  in  LT 
directly,  even  though  LT  is  not  available  according  to  the  definition. 
Not  only  is  a register  saved,  but  also  steps  are  saved. 

(2)  Low  limit  of  DO  loop; 


11= . 


DO.  . . I = II, 


Use  index  register  assigned  to 
I for  II. 


Case  (2)  occurs  when  the  lower  limit  is  an 
expression.  Since  FORTRAN  IV  does  not  allow  expressions  as  DO  limits, 
the  programmer  must  resort  to  the  use  of  a temporary  variable.  The 
system  essentially  undoes  this  wasteful  operation. 


(3)  Upper  limit  of  DO  loop: 
I2=.  . . 


DO. 


I=. 


12, 


Use  register  00  for  12 
if  it  is  available  within 
the  loop. 


The  system  uses  register  00  to  store  the  contents  of 
the  upper  limit  at  the  end  of  the  DO  loop  in  preparation  for  a transfer 
check  back  to  the  beginning  of  the  loop.  As  in  case  (2) , the  system 
essentially  allows  the  upper  limit  to  be  an  expression. 

(c)  Assignment  of  local  variables — local  variables  are 
assigned  to  auxiliary ^registers,  if  available,  in  the  following  order 
(each  is  based  on  decreasing  frequency  of  references) : (1)  variables 
whose  every  reference  is  included  within  a local  domain  and  (2) 
variables  for  which  at  least  one  reference  does  not  appear  in  a local 
domain. 
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A distinction  must  be  made  between  cases  (C-1)  and  (C-2) 
because  the  variable  in  case  (C-1)  is  always  redefined  at  the  beginning 
of  a local  domain.  Hence ^ it  makes  no  difference  if  the  auxiliary 
register  used  to  store  the  variable  value  is  reused  for  another  purpose 
between  two  local  domains  of  the  same  variables.  However,  for 
case  (C-2) , it  is  necessary  to  add  a special  store  instruction  each  time 
an  assignment  is  made,  so  that  a global  value  can  be  maintained  for  the 
varicible . 


It  is  possible  that  a variable  can  be  defined  as  a local 
variable  in  many  local  domains.  Hence,  in  one  domain,  the  variable 
could  be  assigned  register  10  and  in  another  domain,  register  2.  The 
local  variable  register  assignment  is  perhaps  the  most  heavily  used 
portion  of  the  entire  register  assignment  algorithm,  since  the  local 
domains  appear  in  such  frequency  and  are,  in  general,  extremely 
restrictive  in  domain  size.  Cases  (b-1) , (b-2) , and  (b-3)  for  temporary 
variables  apply  also  for  local  variables.  The  register  used  is  no 
longer  availcible  in  the  entire  local  domain. 

The  following  example  should  help  to  illustrate  the 
concept  of  local  variable  register  assignment:  Assume  that  X is 
assigned  registers  1 and  2 in  each  domain,  Y is  assigned  register  2 in 
its  first  domain,  and  Z is  assigned  register  8. 


Code 


X = 1 

E 1 

Y = X+1 

ST  1 

. 

E2 

Z = Y+2 

RE  1 

E4 

ST2 

ST  2 

STORE 

V 

GO  TO  20 

El 

X 

• 

+2 

MARK 

Y = 4 

STORE 

Y 

20 

E2 

ST  LT 

Y = 2+Y 

ST  8 

RECALL 

E 2 

Y 

+8 

XLT 

X = 2 

SEARCH 

STORE 

20 

Y 
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The  variable  X satisfies  case  (c-1)  since  it  is  totally 
defined  within  every  domain  where  it  appears.  Hence,  a register  can  be 
simply  substituted  for  X for  each  reference.  However,  Y appears  as 
case  (c-2) . Its  first  domain  is  local,  but  not  the  remaining 
references.  The  programmer  can  substitute  a register  for  Y in  the 
first  domain,  but  this  register  may  not  be  available  at  statement  20. 
Hence,  one  must  also  save  the  value  for  Y in  its  general  storage 
location . 


(d)  Assignment  of  free  register  variable — the  programmer  can 
specify  that  the  free  registers  be  used  for  variable  assignment  (AUTO 
option,  tcible  I)  . They  are  allocated  in  order  of  decreasing  frequency 
of  references  to  all  variables  not  already  assigned  registers  by  the 
system.  Only  one  variable  can  be  assigned  to  a given  free  register. 

Since  the  compiler  is  designed  to  store  arguments  of 
subroutines  in  basic  registers,  care  must  be  taken  in  s\±>program  calls 
to  save  the  values  of  variables  used  for  permanent  storage  in  nonbasic 
registers  if  they  would  be  destroyed  by  the  called  program.  Hence,  the 
best  feature  of  the  common-expression,  temporary  variable,  and  local 
variable  register  assignment  is  that  the  values  in  the  registers  can  be 
destroyed  when  calls  are  made.  However,  if  free  registers  are  assigned, 
their  contents  may  have  to  be  saved.  The  programmer  does  have  control 
in  specifying  which  registers  need  be  saved  when  the  program  calls  a 
given  subprogrcim  (command  S,  table  IV)  . 


5.1.3  Example;  Use  of  Auxiliary  Registers 

The  sample  program  of  figure  2 illustrates  the  use  of 
auxiliary  registers.  A register  map  is  given  in  figure  3.  The  ARC 
colvunn  gives  the  assigned  register  number  (20  indicates  that  a nonbasic 
register  must  be  assigned) . A maximvim  of  three  work  registers  and  one 
DO  index  register  is  used.  (Actually,  there  are  only  two  registers  after 
B(I)  is  set  up  as  a common  subexpression.  See  A(I)  statement.)  In 
this  example,  all  the  variables  qualify  as  temporary  variables,  and  the 
variables  E,  F,  eind  X are  assigned  registers  00,  14,  and  13, 
respectively.  The  variables  C,  M,  and  K are  strictly  temporary 
variables  that  can  be  assigned  register  15.  The  temporary  variable  N 
satisfies  the  lower  DO  index  condition  and  hence  is  assigned 
register  01.  The  temporary  variable  KK  satisfies  the  upper  DO  index 
condition  and  is  assigned  register  00.  The  variables  A and  B are 
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assigned  nonbasic  registers.  The  variable  J,  although  temporary,  cannot 
be  assigned  an  auxiliary  register  over  its  domain,  since  none  is 
available.  However,  its  domain  can  be  subdivided  into  two  local 
domains,  each  of  which  satisfies  the  strictly  temporary  variable  case  (J 
is  assigned  register  15) . Since,  in  general,  J can  be  assigned  a 
different  register  in  each  domain,  its  ARG  column  has  the  value  20. 


FORTRAN  COMPILER  TEST! 


C PROGRAM  ILLUSTRATES  THE  USE  OF  AUXILIARY  REGISTERS 
DIMENSION  A(l) , B(l) 

E=1 

F=1 

C=1 

X=E+F*C 

M=2 

J=l+M 

K=2*J 

N=K+1 

DO  10  l=N,5 

A(I)=B(I)  + 1./(B(I)*B(D) 

10  CONTINUE 
KK=B ( I ) 

DO  20  1=1 ,KK 
X=X+1 

20  CONTINUE 
J=3 
X=J*J 
STOP 
END 


Figure  2.  Sample  program  listing. 


30 


NAME  TABLE 


NDEX 

NAME 

DIM 

TYPE 

ARC 

1 

A 

1 

0 

20 

2 

B 

1 

0 

20 

3 

E 

0 

0 

0 

A 

F 

0 

0 

]k 

5 

C 

0 

0 

15 

6 

X 

0 

0 

13 

7 

M 

0 

1 

15 

8 

J 

0 

1 

20 

9 

K 

0 

1 

15 

10 

N 

0 

1 

1 

1 1 

1 

0 

1 

20 

12 

KK 

0 

1 

0 

Figure  3.  Register  table. 


It  is  assumed  that  the  loader  will  store  the  pointer  value 
n-1  in  the  register  allocated  to  A if  register  n contains  A(J). 

5.1.4  Program  Execution 

The  relative  Wang  code  generated  is  shovm  in  figure  4 . The 
program  illustrates  the  following  points: 

(a)  Remembering  constants — Steps  0 to  3 store  1 into  E,  F, 
and  C.  However,  since  X has  been  assigned  register  13,  the  system 
changes  ST  15  to  St  13  to  save  a step  and  computes  X in  register  13. 

(b)  Expression  optimization — Steps  4 to  6 translate  X=E+F*C. 
The  parse  string  is  XCF*E+=.  Since  the  system  remembers  that  P is  in 
the  window,  it  does  not  recall  it,  but  simply  multiplies  the  value  by 
the  result  of  X obtained  so  far  (step  4) . 

(c)  Strictly  temporary  variables — Steps  7 to  10  treinslate 
M=2;  J=l+M.  Steps  11  to  13  translate  K=2*J;  N=K+1.  Even  though  M,  J, 
and  K are  using  register  15,  the  register  number  is  changed  to  01  to 
correspond  to  the  first  assignment  for  N.  Step  11  illustrates  the  use. 
of  the  special  "2*"  operator  that  does  a repeated  addition  to  save  a 
step. 
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STEP  CODE  BUTTON  KEY 


STEP  CODE  BUTTON  KEY 


0 

0001 

E 1 

35 

0005 

E 5 

1 

0600 

ST 

00 

36 

0600 

ST 

00 

2 

0614 

ST 

14 

37 

0001 

E 1 

3 

0613 

ST 

13 

38 

0201 

+ 

01 

4 

0413 

X 

13 

39 

0902 

SHIFT 

ALPHA 

5 

0700 

RE 

00 

40 

0806 

SIN 

6 

0213 

+ 

13 

41 

0800 

SEARCH 

7 

0002 

E 2 

42 

DOO 

8 

0001 

ST 

01 

43 

0001 

E 1 

9 

0001 

E 1 

44 

0615 

ST 

15 

10 

0201 

+ 

01 

45 

0801 

RECALL 

11 

0201 

+ 

01 

46 

B 

12 

0001 

E 1 

47 

0215 

+ 

15 

13 

0201 

+ 

01 

48 

1511 

INDIR 

14 

0900  SHIFT 

MARK 

49 

0715 

RE 

15 

15 

DOO 

50 

0912 

SHIFT 

INT 

16 

0615 

ST 

15 

51 

0600 

ST 

00 

17 

0801 

RECALL 

52 

0001 

E 1 

18 

B 

53 

0601 

ST 

01 

19 

0215 

+ 

15 

54 

0900 

SHIFT 

MARK 

20 

1511 

INDIR 

55 

DOl 

21 

0715 

RE 

15 

56 

0001 

E 1 

22 

0600 

ST 

00 

57 

0213 

+ 

13 

23 

0701 

RE 

01 

58 

0001 

E 1 

24 

0615 

ST 

15 

59 

0201 

+ 

01 

25 

0801 

RECALL 

60 

0902 

SHIFT 

ALPHA 

26 

A 

61 

0806 

SIN 

27 

0215 

+ 

15 

62 

0800 

SEARCH 

28 

0700 

RE 

00 

63 

DOl 

29 

0614 

ST 

14 

64 

0003 

E 3 

30  0812 

)(**2 

65 

0615 

ST 

15 

31 

0815 

l/x 

66 

0812 

X**2 

32 

0214 

+ 

14 

67 

0613 

ST 

13 

33 

1511 

INDIR 

68 

0903 

SHIFT 

STOP 

34 

0615 

ST 

15 

69 

0914 

SHIFT 

END 

Figure 

4. 

Wang  code 

listing. 
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(d)  The  IX)  loop  ledael — Steps  14  to  15  define  the  DO  loop  DOO. 
The  system  knows  that  the  index  N is  immediately  available  and  hence  can 
store  it  in  the  first  available  register  used  by  computing  B(I). 

(e)  Common  subexpression — Steps  16  to  22  are  used  to  compute 
the  common  svibexpression  B(I)  and  store  it  in  00. 

(f)  Squaring — Step  30  is  used  to  find  B(I)*B(I). 

(g)  End-of-DO-loop — Steps  35  to  42  perform  the  end-of-DO  loop 
check  (a,  SIN  causes  a two-step  jump  if  the  contents  of  register  00  are 
less  than  the  contents  of  the  display  window) , 

(h)  Upper  DO  index — Steps  43  to  51  tramslate  KK  = B(I)  , 
storing  the  result  in  register  00.  Doing  so  saves  the  updating  of 
register  00  at  the  end  of  the  DO  loop  (see  steps  35  to  36  for  the 
previous  DO  loop) . 

In  many  test  cases,  as  much  as  20  percent  of  the  core  was 
saved  due  to  the  register  optimization. 

5.2  Additions  to  the  Optimizer  Subsystem 

The  second-pass  optimizer  (sect.  3.5)  contains  eight  additional 
passes  through  the  translation  code  in  an  effort  to  reduce  the  number  of 
steps : 


(a)  This  pass  effectively  changes  the  FORTRAN  statement, 
IF(A)N1,N1,N2  to  IF(A.LE.O.)GOTO  N1 

GOTO  N2 

unless  statement  N1  follows.  In  that  case,  IF(A.GT.0.)  GOTO  N2  is  used. 
Also,  the  statement 

IF(A)N1,N2,N2 


is  changed  to 


IF(-A.GT.0. )GOTO  Nl  if  Statement  N2  follows 


or 

IF(-A.LE.O.)GOTO  N2  if  statement  Nl  follows. 
If  neither  follows,  no  change  is  made. 
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(b)  This  pass  changes  the  sequence 


■ Jif+  ■ 

”ch  sign" 

search 

Jif+ 

LI 

search 

search 

to 

L2 

L2 

Mark 

Mark 

-LI 

LI 

IS  ^ 

It  does  this  before  the  pass  that  eliminates  unreferenced  marks,  so  that 
Mark,  LI  is  eliminated  also  if  it  is  not  referenced.  The  code  arises 
from  the  FORTRAN  statement  IF (A),  LI,  L3,  L2  with  statement  LI 
following. 

(c)  This  pass  checks  to  eliminate  trailing  zero  digits  in  a 
number  when  advantageous. 

(1)  If  there  is  only  one  trailing  zero  of  an  integer,  no 
change  is  made. 

(2)  Two  or  more  trailing  zeros  of  an  integer  are  changed 
to  exponential  notation;  for  example,  1000  is  changed  to  1.E3  with  the 
translated  code  El,  a,  f3. 

(3)  Trailing  zeros  of  decimal  numbers  are  eliminated;  for 
example,  12.7300  is  changed  to  12.73. 

(4)  If  the  first  trailing  zero  of  a decimal  number  is 
left  of  the  decimal  point,  the  decimal  point  is  eliminated,  and 
steps  (1)  and  (2)  apply;  for  example,  12500.00  is  changed  to  125E2, 
with  the  a shift  used  to  multiply  by  100. 

(5)  All  the  trailing  zeros  of  a number  with  an  exponent 
to  the  right  of  a decimal  point  are  eliminated;  for  example,  12.50E23  is 
changed  to  12.5E23.  This  change  does  not  result  in  an  optimum 
translation,  since  the  decimal  point  can  be  eliminated  by  use  of  125E22. 
It  is  eliminated  in  pass  (d) . 

(6)  If  the  first  trailing  zero  of  a number  with  an 
exponent  is  to  the  left  of  a decimal  point  or  just  to  the  right  of  it, 
the  decimal  point  is  eliminated,  and  the  exponent  is  adjusted;  for 
example,  -120.00E25  is  changed  to  -12E26,  and  32.0E-15  is  changed  to 
32E-15. 

(7)  The  code  EO,  ST  R,  is  changed  to  T R. 
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The  net  effect  of  this  pass  is  to  make  sure  that  no  number 
is  written  with  a trailing  zero  (in  the  mantissa,  if  in  exponential 
notation)  unless  it  is  an  integer  multiple  of  10,  but  not  of  10  for 
n > 2. 


(d)  This  pass  insures  that  identical  numbers  have  identical 
representations . 

(1)  Integers  are  left  untouched. 

(2)  Decimal  points  are  introduced,  if  an  exponent  can  be 
eliminated;  for  example,  12E-1  is  best  written  1.2. 

(3)  An  exponent  of  1 is  eliminated  for  a nondecimal 
mantissa;  for  example,  12E1  is  changed  to  120. 

(4)  A SETEXP  notation  is  changed  to  an  a shift  by 
introduction  of  a decimal  point;  for  example,  12E-16  is  best  written 
1.2E-15.  The  exponent  field  can  now  be  translated  in  two  steps  (a,  F15) 
instead  of  four,  with  a net  saving  of  one  step. 

(5)  All  leading  zeros  of  decimal  numbers  are  eliminated; 
for  example,  0.03  is  represented  as  3E-2  and  0.0001  as  lE-4. 

(6)  For  a decimal  point  with  an  exponent,  the  exponent  is 
eliminated,  if  possible,  by  shifting  the  decimal  point  to  left;  for 
excimple,  12.3E-2  is  changed  to  .123. 

(7)  Otherwise,  a check  is  made  to  see  if  a SETEXP 
notation  can  be  changed  to  an  o shift;  for  example  12.3E-16  is  best 
written  1.23E-15  as  in  item  (4) . 

(8)  Otherwise,  a check  is  made  to  see  if  the  exponent  can 
be  eliminated  moving  the  decimal  point  to  the  right;  for  excimple, 
12.36E2  is  changed  to  1236,  and  0.013E2  is  changed  to  1.3. 

(9)  Otherwise,  a check  is  made  to  see  if  the  decimal 
point  can  be  eliminated  by  moving  it  to  the  right;  for  example,  12.3E5 
is  changed  to  123E4  (but  not  12.3E-15  to  123E-16) . See  item  (4)  about  A 
SETEXP  notation. 

The  net  effect  of  passes  (c)  and  (d)  is  to  obtain  a unique 
translated  representation  for  all  numbers. 

(e)  This  pass  eliminates  the  necessity  of  entering  a number 
when  it  is  determined  that  it  is  in  the  window.  Thus,  the  FORTRAN  code 

A(l)=1.2 

X=1.2 

Y=1.2 
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is  translated 

generated  [A(D]  in  LT 
E 1 
E . 

E 2 
INDIR 
ST  LT 
STORE 
STORE  ^ 

(f)  This  pass  further  attempts  to  optimize  steps  involving 

arrays . 

(1)  If  the  value  of  a needed  array  is  in  the  window,  it 
is  not  generated  again.  For  example, 

A(I)  = . . . 

T = A(I)+. 


is  obvious. 


(2)  If  the  index  of  an  array  is  determined  to  be  in  a 
register,  it  is  not  generated  again.  For  example, 

A(I)=  .... 

X = Y 

T = A(I)+..  . . . 

The  determination  is  made  that  £[A(I)]  is  still  in  LT,  so  the  second 
A (I)  is  translated  INDIR,  RE  LT.  Another  example  is 

T = B(I)*A(I) 

B(I)=X 

A(I)=Y 

where  T is  a basic  register.  The  computation  is  performed  in  that 
register,  leaving  i[A(I)l  and  i[B(l)]  undestroyed. 

(3)  Item  (2)  applies  if  the  index  differs  only  by  a 
constcint.  A code  to  subtract  or  add  the  constant  is  supplied  as  needed. 
This  pass  was  motivated  by  frequent  use  of  em  exchange  code  in  problems 
involving  arrays.  Thus,  the  code 


T=A(I) 

A(I)=B(I) 

B(I)=T 
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is  now  translated  with  the  locations  of  A(I)  and  B(I)  generated  only 
once. 


(g)  This  pass  tries  to  avoid  an  extraneous  code  involving  the 
use  of  register  00. 

(1)  It  eliminates  T 00  when  not  needed.  For  excimple, 

IF (A) 5,5,10 
5 IF(B)15,15,20 

15  .... 

produces  a code  involving  T 00  twice.  The  second  T 00  is  eliminated. 

(2)  It  eliminates  the  code  RECALL  N,  ST  00  when  not 
needed,  such  as  occurs  in  the  translation 

D05I=1,N 

D05J=I,N 

(3)  It  replaces  the  code  RE  00,  . . .,  T 00,  by  T 00, 

. . . provided  that  the  user  has  not  specified  a COMPILE  R for  00  or 
that  the  meaning  has  not  changed.  This  item  is  for  situations  in  which 
00  is  used  as  a temporary  varieible  prior  to  a zero  being  needed  in  00 
for  any  reason. 

(h)  A final  check  is  made  to  eliminate  all  recalls  of 

nondimensional  variables  that  are  in  the  window.  If  an  RE  R is  directly 
preceded  by  PRINT,  format,  or  INDIR,  ST  R,  or  any  jump,  SEARCH,  Ll,  or 
combinations  of  these,  and  the  next  previous  operation  is  a pR,  then  the 
RE  R is  eliminated.  Similarly,  if  the  code  is  STORE,  X or  RECALL,  Y 
followed  by  the  above  nondestructive  steps,  followed  by  RECALL,  X,  the 
last  two  steps  are  eliminated.  After  any  jump,  a STOP  or  RETURN  is 
allowed  in  place  of  SEARCH. 

6.  EXAMPLE 

This  lengthy  example  illustrates  the  use  and  effectiveness  of  this 
compiler.  The  example  consists  of  a main  program  that  calls  the 
subroutine  POLRT,  teJcen  without  change  from  the  360  Scientific 

Subroutine  Package.^  The  POLRT  is  used  to  find  the  complex  roots  of 


^Syst&n/360  Scientific  Subroutine  Package,  Version  III,  Programmer's 
Manual,  Program  Number  360A-CM-03X,  ed.  GH20-0205-4 , IBM  Corporation, 
White  Plains,  NY  (August  1970),  181-183. 
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real  polynomials;  the  calling  program  must  supply  the  degree  M and  the 
M+1  coefficients  XCOF  of  the  polynomial.  A simple  main  program  was 
therefore  written  to  read  these  inputs,  call  POLRT,  and  write  the 
output.  Figure  5 shows  the  program,  with  comments  to  give  motivations 
for  each  of  the  COMPILE  statements  that  affect  the  translation.  The 
MAIN  was  made  a subroutine  to  allow  the  arguments  to  be  aligned 
identically  to  POLRT.  This  is  a common  maneuver  to  save  steps  when 
subroutines  are  called,  but  it  does  require  the  COMPILE  S statement. 
The  dimensions  of  arrays  were  guessed  at  to  try  to  allow  large-degree 
polynomials  in  a fully  expanded  Wang  520. 

In  the  *FTC  card  for  MAIN,  no  options  were  specified;  hence,  both 
the  source  code  and  the  translation  output  (name  table,  translation,  and 
reference  table)  are  printed.  Figures  6 to  8 show  the  translation 
output.  These  have  been  previously  described,!  except  that  the 
reference  table  (fig.  8)  now  contains  three  additional  lines:  The  top 
work  register  is  09  (registers  10  to  15  contain  the  arguments) ; only 
single  DO  loops  are  involved  (because  of  reading  and  writing  arrays) ; 
and  one  is  the  maximum  number  of  work  registers  needed  to  translate  any 
statement . 

In  the  *FTC  card  preceding  POLRT,  the  AUTO  option  was  specified. 
Since  POLRT  calls  no  other  routines  and  all  of  its  nonargument  variables 
are  initialized  at  every  entry,  many  steps  are  saved  by  allowing  the 
compiler  to  assign  basic  registers  for  the  most  frequently  used 
varicible.  If  some  of  its  nonargument  variables  were  not  initialized 
every  time,  or  if  POLRT  had  been  called  in  a loop,  then  some  register 
conflicts  could  arise.  Then  AUTO  (nn-mm)  would  have  to  be  used  by  the 
programmer  to  specify  a more  limited  range.  It  is  always  up  to  the 
programmer  to  insure  that  no  conflicts  arise  when  many  routines  are 
compiled  together. 

Figure  9 contains  the  source  listing  of  POLRT  (the  original  comments 
have  been  deleted) ; figure  10,  the  name  table;  figure  11,  the 
translation;  and  figure  12,  the  reference  table.  In  figure  12,  the  work 
registers  are  09,  08,  and  07,  and  register  01  is  used  for  the  single  DO 
loops.  If  AUTO  were  not  specified,  registers  02-06  would  never  be  used. 
The  name  table  in  figure  10  assigns  N,  X,  Y,  SUMSQ,  and  IFIT  to  these 
registers.  In  addition,  the  compiler  senses  that  TEMP  can  be  assigned 
to  one  of  the  work  registers  (07) , since  it  is  a temporary  variable  in  a 
single  domain  that  does  not  use  register  07.  Similarly,  the  variable  DX 
was  assigned  register  01,  because  it  is  defined  and  used  in  a single 
domain  outside  a DO  loop,  and  the  variables  XTZ,  L,  and  ITEMP  were  all 
assigned  register  00.  (When  L is  used  as  an  index,  it  is  assigned 


^H.  Bloom  and  A.  Hausner,  FORTRAN  IV  Compiler  for  the  Wang  520 
Calculator,  Harry  Diamond  Laboratories  TM-73-15  (July  1973), 
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FORTRAN  COMPILER  MAIN 


SUBROUTINE  MA I N ( XCOF , CGF , M, ROOTR , ROOTI , I ER) 

DIMENSION  XCOF(25)tCOF(25),ROOTR(2A),ROOTK24J 
REAOn  M 

COMPILE  W M{  I503r 1,31 

C THIS  W CARD  PRINTS  M ON  THE  READ  IN  FORMAT  1500  WITH  SPACES  BEFORE' 

C AND  AFTER  THE  PRINT. 

N = M + l 

READO  (XCOF(  n,  I = l,N) 

COMPILE  W XCOFI05 10, 1,0) 

C THIS  W CARD  PRINTS.  THE  COEFFICIENTS  XCOF  ON  THE  READ  IN  FORMAT  0510. 

CALL  POLRT(XCOF,COF,M,ROOTR,ROOTI ,IER) 

COMPILE  S POLRTlE,E,E,E,E,E  ),0-0 

C THIS  S CARO  INDICA/ES  THAT  NO  ARGUMENTS  OF  MAIN  NEED  BE  SAVED  OR 
C RESTORED,  AND  THAT  NO  ARGUMENTS  OF  POLRT  NEED  BE  INPUTTED  OR  OUT- 
C PUTTED.  THE  REASON  MAIN  WAS  MADE  A SUBROUTINE  WAS  TO  MATCH  ARGUMENTS 
C FOR  THIS  PURPOSE. 

IF(  lEP.NE.O)  GOTO  5 

WRITE!)  IROOTR! I ) ,ROOTI( n , 1=1, M) 

COMPILE  W ROOTR ( 0010, 0, I ) 

COMPILE  W ROOTI ( i 110,0,0 ) 

C THESE  W CARDS  PRINT  THE  REAL  AND  IMAGINARY  ROOTS  IN  FORMATS  0010  AND 
C 1110,  SEPARATING  EACH  PAIR  OF  ROOTS  BY  SPACES. 

STOP 

5 WRITE!  ) lEP 
COMPILE  W lER  !0700,0, 1) 

C THIS  W CARO  PRINTS  THE  VALUE  OF  lER  IN  THE  FORMAT  0700  AFTER  A SPACE. 
STOP 
END. 


Figure  5.  Program  illustrating  use  of  COMPILE  statements 


tablf 


INDEX  NA«E 


DIM  TYPE  ARG 


1 MAIN 

2 XCPE 

3 COF 

4 M 

5 ROOTR 

6 - ROOT  I 

7 lER 

8 N 

9 I 

13  POLRT 


0 I 50 

25  0 15 

25  0 14 

3 I 13 

24  0 12 

24  0 11 

0 I 10 

0 I 20 

0 I 20 

3.0  50 


Figure  6.  Ncune  table  for  program  MAIN 


register  01.)  These  temporary  variables  would  have  been  assigned 
auxiliary  registers  even  if  AUTO  had  not  been  specified.  The  net  effect 
of  AUTO  was  just  the  assignment  of  variables  in  registers  02-06. 

When  the  reader  examines  the  translation  of  POLRT  in  figure  11,  he 
must  keep  in  mind  some  of  the  new  optimization  features  automatically 
used.  For  example,  N+1  is  a common  subexpression  in  the  line  below 
statement  number  35,  NXX=N+1  and  two  lines  following  KJ1=N+1.  Register 
00  is  used  for  this  subexpression,  computed  in  steps  56  to  58,  and  later 
recalled  in  step  64.  Also,  KJl  is  stored  in  register  07  (step  65)  for 
use  in  the  DO  loop  following  (steps  72  and  88) . Register  07  was 
available  for  that  local  domain  and  one  step  was  saved  by  this  action. 
(If  KJl  had  been  a basic  register,  one  step  would  have  been  wasted. 
This  feature  is  implemented  before  register  assignments  and  does 
sometimes  lead  to  wasted  steps . ) 

Following  the  POLRT  source  deck  was  an  *LDR  card  and  an  *END  card 
with  no  options  listed.  The  storage  maps  (fig.  13)  and  progreun  listing 
(fig.  14)  are  printed  regardless  of  any  options  that  affect  only  entry 
point  codes,  mark  numbers,  and  registers  assigned.  In  figure  13,  the 
arrays  were  assigned  registers  16  to  113  as  listed  under  LIMITS  in  the 
storage  map  for  MAIN.  (These  are  not  repeated  in  the  storage  map  for 
POLRT,  because  the  dimensions  there  were  given  as  one.  All  varicible 
dimensioned  arrays  should  have  dimension  one.)  The  other  columns  are 
almost  self  explanatory.  The  INDEX,  NAME,  DIM,  TYPE,  and  ARG  columns 
are  repeated  from  the  name  table  (fig.  6)  only  for  variables  that  must 
be  assigned  space.  Columns  for  COM  and  EQU  are  for  variedsles  in  common 
or  equivalence.  The  LOC  column  contains  the  location  of  the  variable 
or,  in  the  case  of  an  array,  the  location  of  the  pointer.  These 
locations  can  be  used  to  recall  these  variables  with  the  indirect  code. 
To  recall  varieibles  with  the  direct  code,  the  REG  column  lists  the 
direct  register  name.  The  LIMITS  colxomn  is  applicable  only  to  arrays. 
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EP 

CODE 

BUTTON 

KEY 

STEP 

CODE 

BUTTON 

KEY 

0 

0900 

SHIFT 

MARK 

45 

0804 

J IF  0 

1 

MAIN 

46 

0800 

SEARCH 

2 

0802 

PR  INT 

47 

1 

5 

3 

0015 

CLRDSP 

48 

0001 

E 1 

4 

0004 

E 4 

49 

0601 

ST 

01 

5 

0609 

ST 

09 

50 

0900 

SHIFT 

• MARK  • 

6 

0903 

SHIFT 

.STOP 

51 

L ABOl 

7 

0613 

ST 

13 

52 

0802 

PR  INT 

S 

0302 

/ 

PR  INT 

53 

0015 

CLRDSP 

9 

1500 

SP-RE 

00 

54 

0701 

RE 

01 

10 

0802 

PRINT 

55 

0609 

sr 

39 

11 

0015 

CLRDSP 

56 

0712 

RE 

12 

12 

0609 

ST 

09 

57 

0209 

’ *■ 

09 

13 

0001 

E 1 

58 

1511 

I NDI  R 

14 

0209 

+ 

09 

59 

0709 

RE 

09 

15 

0901 

SHIFT 

STORE 

60 

0802 

PRINT 

16 

N 

'61 

0010 

E . 

17 

0001 

E 1 

62 

0701 

RE 

01 

18 

0601 

ST 

or 

63 

0609 

ST 

09 

19 

0900 

SHIFT 

MARK 

64 

0711 

RE 

11 

20 

LABOO 

65 

0209 

+ 

09 

21 

0609 

ST 

09  ' 

66 

1511 

INDIR 

22 

0902 

SHIFT 

ALPHA 

67 

0709 

RE 

09 

23 

1103 

SHFT  F 

03 

68 

0802 

PRINT 

24 

0600 

ST 

00 

69 

1110 

SHFT  F 

10 

25 

0715 

RE 

15 

70 

0713 

RE 

13 

26 

0209 

<■ 

09 

71 

3600 

ST 

33 

27 

0002 

E 2 

72 

0001 

E 1 

28 

0200 

00 

73 

0 201 

+ 

01 

29 

090  3 

SHIFT 

STOP 

74 

0902 

SHIFT 

ALPHA 

30 

1511 

IND  IR 

75 

0806 

SIN 

31 

0609 

ST 

09 

76 

3800 

SEARCH 

32 

0802 

PRINT 

77 

LABOl 

33 

0510 

/ 

10 

78 

0903 

SHIFT 

STOP 

34 

0801 

RECALL 

79 

0930 

SHIFT 

MARK 

35 

N 

80 

5 

36 

0600 

ST 

00 

81 

38  02 

PRINT 

37 

000  1 

E ' 1 

8 2 

0015 

CLRDSP 

38 

0201 

+ 

01 

83 

0710 

RE 

10 

39 

0902 

SHIFT 

ALPHA 

84 

3832 

PRINT 

40 

0806 

Sin 

85 

0700 

RE 

00 

41 

0800 

SEARCH 

86 

090  3 

SHIFT 

STOP 

42 

L ABOO 

37 

0914 

SHIFT 

END 

43 

POLPT 

44 

0710 

RE 

10 

Figxire  7.  Wang  relative  code 

for  program 

MAIN. 
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REFERENCE ' tAB?^  M 


MAIN 

N 

POLRT 

5 ^ 

LABbb 


. • ' ^ --V?  * ' 

’vc-*; 

I ' 

••  Rl 

j'J.  Ji,' 


Tv  ■•.'•r,  ,.,v  •, 


'.1 

16 

*3 


35 


f ui RE^;;-- i; ,,  4.  £“.r:?3  - 

® ’V- yV #;■'.•  . V •.  • •..•  ’•  f 4-f 

■ AD'Ah'  •'  '^0 


•*0-» 


THE  rrilsi^  AviTtAmE’  is  9 

NAX  Ob  LOOP  NEST  is  I ^ 

HAX  REG.liANGE  JS  I ' '!  « • 

Figure  8.  Reference  table  for  ppogreim  MAIN, 


FORTRAN  COMPILER  POLRT 


AUTO 


^ SUBROUTINE  POLRTIXCOFfCOF.M.ROOTRiROO'TIilERI  f ^ ■ 

0 IHENS lON^XCOF II » , COFI I » , ROOTRI  ij  ,R00tl < H # > > ' 

DOUBLE  PRECISION  X6t VOt'X, V, XPRVTPR“fUX,UVt.Vf YTrxT.Ut^TE 
I DX.DV,* TEMP, ALPHA 
iFlTiQf  V 

. iER*o.  ' 

|F(XCOF(Nflt)10,25,lO 
“ 10  T FIN  I 1 5;  is;,  32 
15  IER*l 
20  RETURN 
i?  iERi'v:f^'"-i:/.-. 

‘ -6Q TO'  20,.  /,  f.; 

I?  GO  TO  20 
;:jii2  IFIN-361  35,35,30 
Ti5  NX»N 


‘NXX»N«-1 


:>$;mi v,, 

""  v,bO  vW  , k3i  ‘ • 
Mf*KJi-L>i  . 

To  COFiMtI'XCOFItl 
T$;'x'0*rpo  5001  oT  ^ ; 
»vl"  VO*O.0lOOOlOl 
V^  -IN^O  . V. 

p;  x-xo:;^o;  ■*:' 

L#.xo*-to.otyo.' 

5St^l¥0?rl,0.0#X,'  ' 


Figure  9.  Program  listing  for  POLRT. 
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X=XO'  “ 

Y*YO 
IN>IN4-l 
GO  TO  59 

55  IFIT=l 
XPR*X 
YPR*Y 

59  ICT»0 

60  UX=0.0 
UY=0.0 
V =0.0 
YT*0.0 
XT=l.O 
U=C0F(N»14 
IFtU»  65.130,65 

65  00  70 
L =N-I«-l 
TEMP=COF(L  » 
XT2=X*XT-Yil'yT 
YT2=X*YT+Y*XT 
U«U«-TEMP*XT2 
V=V»TEMP*YT2 
FI  = I 

UX=UXtFI*XT*TEMP 

UY=UY-FI*YT*TEMP 

XT=XT2  . 

70  YT=YT2 

S 0 MS  Q=U  X^'O  X f UV*  Uf 
IFtSUMSO)  75,110,75 

75  OX=(V*UY-U*UX)/SUMSO 


DY=- IU*UY*V*UX 1/SUMSQ 
Y=Y>OY 

78  IFlDABSlOYl+O^kBSlOXl-L. 00-05)  100,80,80 
80  ICT=lCT*l 

lFCICT-500)  60,85,85 
■"85  I F n F i'T ) 100^90  , 100 
90  lFllN-5)  50,95,95 
95  IER=3  /. 

GO  TO  20 

100  DO  105  L*l,NXX 
MT=KJl-L*l 
"TEMP=xCOFfMTr' 

XCOF(MT)=COFlL) 

105  COFtL)«TEMP  , 

ITEMP=N 

N=NX. 

NX=ITEMP 

■ • • I F (I F I T ) i 2 Of  5 5; 1 20 
no  IFIIFIT)  11^,50,115' 

115  X=XPR^. 

Y=YPR 

120  !FIT=0 

122  IFIDABS FY1-l.0D-4*DABS(X»l  135,125,125 
Figure  9.  Program  listing  for  POLRT  (cont'd) . 
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125  ALPHA=X+X 

SUMSO=X*X+Y«Y 

N=N-2  _ - - 

GO  TO  lAO  , 

130  X=0,0 
NX=N<-1 
NXX=NXX-l 
135  Y=0.0 

SUMSQ=0.0 
ALPHA=X 
N = M-l 

lAO  COF(  2)  = C0F  ( 2)  ♦ALPHA=<'CUFl  I » 

145  DO  150  L=2,N 

150  COF(L  ) = COF(L  + l M-ALPHA*COF(L  l-SUMSO*COF(L-U 
155  R00TI(M2)  = Y 
ROOTR (N2)=X 
N2=N2+l 

IF(SJMSO)  160,  165,160 
160  Y=-Y 

SUMS0=0.0 
GO  TO  155 

165  IF(N)  20,20,45 
END 


Figure  9.  Program  listing  for  POLRT  (cont'd) . 


Listings  under  LABEL  and  MARK  indicate  the  mark  code  given  for  each 
label.  In  figure  13  MAIN  used  marks  0000  to  0002,  and  POLRT  used 
marks  0003  to  0107. 

No  conflicting  marks  or  registers  exist,  despite  that  register  114 
was  assigned  to  the  variable  N in  MAIN  and  the  variable  XO  in  POLRT. 
This  is  part  of  an  optimization  scheme  that  equivalences  variables  in 
several  programs.  In  effect,  the  variables  are  placed  in  common.  To  be 
a candidate  for  such  treatment,  the  domain  of  a variable  must  not 
contain  any  function  or  subroutine  calls,  and  the  variable  must  not  be 
in  an  EQUIVALENCE  or  COMMON  statement.  Such  variables  are  considered 
local  to  the  problem  and  can  be  destroyed  after  leaving  the  program. 

The  final  listing  is  shown  in  figure  14.  Array  pointers  are  first 
set  up  (steps  02  to  13).  Next,  the  routines  MAIN  and  POLRT  are  listed 
with  the  assigned  registers,  labels,  and  entry  points  replacing  the 
symbolic  ones.  All  statement  references  in  figxires  8 and  12  are  so 
replaced.  Finally,  the  large  numbers  .5001010-02  and  .10001010-01 
needed  in  steps  99  and  103  of  POLRT  (fig.  11,  12)  are  generated  as  data 
in  the  program  code  in  steps  720  to  735.  Because  any  steps  converted  to 


44 


NAME  TABLE 


INDEX 

• NAME 

DIM.  TYPE 

ARC 

a 

POLRT 

0 

0 

50 

2 

XCOF 

1 

0 

15 

3 

COP 

1 

0 

14 

4-'  ■■ 

-A., 

0 

1 

13 

5 

vROOTR 

r 
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0 

0 

20 

15 
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0 

0 
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27 
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28 
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31 
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1 
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33 
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.0 

Figvire  10.  Name 
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POLRT. 
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STEP  CnOE  •■^UTTCN  KEY  STEP  CODE  nUTTON  KEY 
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PE 
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SHIFT 
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0100 
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0708 
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ST 
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09 

33 
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07 
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39 

*0600 

ST 
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35 
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SHIFT 

MARK 
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3001 
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36 

30 

91 
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01 

37 
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E 2 

92 

Q902 

SHIFT 

ALPHA 

3*8 

0610 

ST 
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93 

0806 

SIN 
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0800 

SEAPCH 

94 

0300 

SEARCH 

40 

20 

95 

DOO 

41 

090C 

SHIFT 

MARK 

96 

0900 

SHIFT 

MARK 

42 

32 

97 

45 

43 

0702 

RE 

0? 

99 

C901 

RECALL 

44 

0609 

ST 

09 

99 

.5001010000000 

45 

0100 

T 

00 

100 

0901 

SHIFT 

STORE 

46 

000  3 

F 3 

lOl 

XO 

47 

0006 

E 6 

102 

0801 
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48 

0309 

- 

09 

103 
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49 

0902 

SHIFT 

ALPHA 
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0901 

SHIFT 
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50 

0305 
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54 
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STORf: 
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SHIFT 
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STEP  Ct'Dh  fjJTTON  KEY  STEP  CODE  BUTTON  KEY 


no 

50 

165 

UY 

1 1 1 

C801 

RECALL 

166 

0901 

SHIP  T 

STORE 

m 

xn 

167 

V 

113 

0608  ST 

03 

163 

0901 

SHIFT 

STORE 

114 

0603  ST 

03 

169 

YT 

115 

0001 

E 1 

1 70 

0001 

E 1 

l 16 

0000 

E 0 

171 

0901 

SHIFT 

STORE 

117 

0012 
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172 

XT 

1 IS 
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00 

1 73 

0702 

RE 

02 

119 

0607  ST 

07 

1 74 

0609 

ST 

09 

m 

0801 

RECALL 

1 75 

0001 

E 1 

121 

YQ 

176 

0209 

f 

09 

122 

0407  < 

07 

1 77 

0714 

RE 

14 

123 

0901  SHIFT 

STORE 

178 

0209 

♦ 

09 

124 

xn 

179 

1511 

INOIR 

125 

0700  RF 

00 

1 90 

0709 

RE 

09 

126 

0601  ST 

01 

ISl 

0901 

SHIFT 
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127 

0703  RE 

03 
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U 

128 

0401  < 

01 

133 

0904 

SHIFT 
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129 

0901  SHIFT 
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134 
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SEARCH 

130 
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130 
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07 
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MARK 
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ST 
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138 

0001 

E 1 

193 
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09 
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140 
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0200 
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143 
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MARK 
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220 

0801 
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XT 

221 
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2 22 
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Figure  11.  Wang  code  for  POLRT  (cont'd). 
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STEP 
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Figvire  11.  Wang  code  for  POLRT  (cont'd). 
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448  50 
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455  0900  SHIFT 

456  120 
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459  0913  SHIFT 
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462  0703  PF 

463  0913  SHIFT 

464  0902  SHIFT 
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466  0309 

467  0012 

468  0902  SHIFT 

469  0805 

470  0800 

471  135 

47?  0703  FE 

473  0609  ST 

474  0209  f 

475  0901  SHIFT 

476  ALPHA 

477  0703  RF 

478  0812 

479  0605  ST 
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Figure  13.  Storage  map  for  MAIN  and  POLRT  (cont*d) 
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0710 

RE 

n 

429 

0705 

RE 

05 

375 

0801 

RECALL 

430 

0509 

/ 

09 

376 

0712 

PE 

12 

431 

0901 

SHIFT 

STORE 

377 

0901 

SHIFT 

STORE 

432 

0713 

RE 

13 

378 

0709 

RE 

09 

433 

0204 

♦ 

04 

379 

07C2 

RE 

02 

434 

0801 

RECALL 

380 

0600 

ST 

00 

435 

0713 

RE 

13 

381 

0001 

E 1 

436 

0913 

SHIFT 

ABS 

332 

0201 

f 

01 

437 

0609 

ST 

09 

303 

0902 

SHIFT 

ALPHA 

438 

0701 

RE 

01 

384 

0806 

SIN 

439 

0913 

SHIFT 

ABS 

Figure  14.  Final  Wang  code  listing  (cont'd). 
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STFP 

CODfr 

OUT  TON 

KEY 

STFP 

CODE 

BUTTON 

KEY 

7209 

f 

09 

495 

03C0 

_ 

03 

441 

7100 

T 

OC 

496 

0001 

E 1 

44? 

0001 

F 1 

49  7 

0200 

f 

00 

443 

090  2 

SHIFT 

ALPHA 

49P 

0609 

ST 

39 

444 

1105 

SHFl  r 

05 

499 

0715 

RE 

15 

445 

:J:>09 

- 

09 

500 

0?C9 

•*' 

09 

446 

0012 

CHS 

501 

1511 

INOI  R 

44  7 

093? 

SHIFT 

Al  PHA 

50  2 

0709 

RF 

09 

448 

0305 

J IF  4- 

503 

060  7 

ST 

07 

449 

0800 

SFACCh 

504 

0700 

RF 

00 

450 

7100 

T 

00 

505 

0609 

ST 

09 

451 

0901 

P FCALL 

506 

0715 

RE 

1 5 

45? 

0304 

SP 

04 

507 

0209 

09 

453 

0 60  9 

ST 

09 

508 

0701 

RF 

31 

4 54 

0001 

r 1 

509 

0608 

ST 

38 

455 

0209 

f 

09 

510 

0714 

RE 

14 

456 

OPOl 

SHITT 

STC-RF 

511 

0208 

f 

08 

457 

0804 

SP 

04 

512 

1511 

INDIP 

458 

0005 

F 5 

51  3 

3 7 03 

RF 

38 

459 

0902 

SHIFT 

Al  PHA 

514 

151  1 

INOIR 

460 

1002 

r ( X ) 

C2 

515 

3609 

ST 

09 

46  1 

0309 

- 

09 

516 

7 707 

RF 

37 

462 

0012 

ChS 

51  7 

1511 

INDIR 

463 

790  2 

SHIFT 

ALPHA 

518 

0603 

ST 

08 

464 

0 305 

J IF  f 

519 

0301 

RECALL 

465 

0300 

SEARCH 

520 

38  0 0 

SP 

33 

466 

0012 

ALL  UP 

1? 

521 

0600 

ST 

00 

467 

0706 

8E 

Ob 

522 

0001 

F 1 

469 

0004 

J IPO 

52  i 

0201 

+ 

01 

46  9 

0800 

SF  ARCH 

52  4 

0902 

SHIFT 

ALPHA 

47J 

JlOO 

T 

00 

5?5 

0806 

SIN 

471 

0301 

recall 

5?t. 

0300 

SEARCH 

472 

030  3 

SP 

03 

52  7 

01  01 

T 

31 

47  3 

0679 

ST 

7P 

523 

770? 

RE 

32 

474 

0005 

F 5 

529 

0600 

ST 

33 

475 

0309 

- 

09 

530 

0301 

RECALL 

476 

0012 

CEiS 

531 

0715 

RF 

15 

477 

0902 

SHIFT 

ALPHA 

532 

0602 

ST 

7? 

478 

0 805 

J IF  f 

533 

0700 

RF 

00 

4 !9 

0800 

search 

534 

0901 

SHIP  T 

STORE 

480 

0 009 

All  UP 

09 

535 

>0715 

RE 

1 5 

481 

7003 

F 3 

53  6 

0706 

RE 

06 

482 

06  10 

ST 

10 

53  7 

0904 

Shift 

J NE  0 

.483 

0800 

SEARCH 

538 

0800 

SEARCH 

434 

0005 

ALl  UP 

05 

537 

001  1 

ALL  UP 

1 1 

405 

0900 

SHIFT 

MARK 

54  0 

0300 

SEARCH 

486 

0100 

T 

00 

5^l' 

no2 

T 

32 

487 

700  1 

r i 

54  2 

0900 

SHIFT 

MARK. 

438 

0601 

ST 

01 

543 

0015 

ALL  UP 

15 

439 

09  n 

SHIFT 

mark 

54  4 

0706 

RE 

36 

49  0 

01  01 

T 

01 

54  5 

0904 

SHIFT 

J IMF  3 

491 

0301 

PlCALL 

54o 

0300 

SEARCH 

47? 

0807 

SP 

02 

54  7 

0009 

ALL  .'JP 

09 

4^3 

0600 

ST 

70 

543 

0801 

RECALL 

494 

0 70  1 

RF 

01 

54  9 

0704 

RE 

04 

Figure  14.  Final  Wang  code  listing  (cont'd). 
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s rFP 

croE 

BUTTHN 

KEY 

STEP 

CODE 

BUTTON 

KEY 

5b0 

0603 

ST 

03 

605 

0104 

T 

04 

5t>l 

0801 

RECALL 

606 

0105 

T 

05 

552 

0705 

RE 

05 

607 

0703 

RE 

03 

553 

0604 

ST 

04 

608 

0901 

SHIFT 

STORE 

55^ 

0900 

SHIET 

MARK 

609 

0714 

RE 

14 

S5S 

0102 

T 

02 

610 

0001 

E 1 

556 

0106 

T 

06 

611 

0302 

- 

02 

557 

0704 

RE 

04 

612 

0900 

SHIFT 

MARK 

558 

0913 

SHIf-T 

AfiS 

613 

0104 

T 

04 

55^ 

0609 

ST 

09 

614 

0002 

E 2 

560 

0100 

T 

00 

615 

0609 

ST 

09 

561 

0705 

RE 

03 

616 

0714 

RE 

14 

562 

09  13 

SHIFT 

ABS 

617 

0209 

♦ 

09 

563 

0902 

SHIFT 

AL  PHA 

6 10 

0001 

E 1 

56^ 

1 104 

SFFT  F 

04 

619 

0608 

ST 

08 

565 

0309 

- 

09 

620 

0714 

RE 

14 

566 

0012 

CHS 

621 

0208 

♦ 

08 

567 

0902 

SHIFT 

ALPHA 

622 

1511 

INDIR 

568 

0805 

J IF  f 

623 

0708 

RE 

08 

569 

3900 

SEARCH 

624 

0608 

ST 

08 

570 

0 10  3 

T 

03 

625 

0801 

RECALL 

571 

C703 

RE 

03 

626 

0714 

RE 

14 

572 

0609 

ST 

09 

627 

0408 

X 

08 

573 

0209 

f 

09 

623 

1511 

INDIR 

574 

090  1 

SHIFT 

STORE 

629 

0209 

f 

09 

575 

0714 

RE 

14 

630 

0002 

E 2 

576 

0703 

RE 

03 

631 

0601 

ST 

01 

577 

08  12 

X**2 

632 

0900 

SHIFT 

(HARK 

578 

0605 

ST 

05 

633 

0105 

T 

05 

579 

0704 

RE 

04 

634 

0701 

RE 

01 

580 

0812 

635 

0609 

ST 

09 

58  1 

0205 

f 

05 

636 

0001 

E 1 

592 

000? 

E 2 

637 

0209 

09 

583 

0302 

- 

02 

638 

0714 

RE 

14 

584 

0900 

SEARCH 

6 39 

0209 

f 

09 

585 

01  04 

T 

04 

640 

0701 

RE 

01 

586 

0900 

SHIFT 

MARK 

641 

0608 

ST 

08 

58  7 

OOli 

ALL  UP 

13 

642 

0714 

RE 

14 

598 

0103 

T 

05 

643 

0208 

08 

C801 

RECALL 

644 

1511 

INDIR 

5Q0 

C7L5 

PE 

15 

645 

0708 

RE 

06 

591 

0609 

ST 

09 

646 

0600 

ST 

08 

592 

OCOl 

E 1 

647 

0801 

RECALL 

593 

0309 

- 

09 

648 

0714 

RE 

14 

594 

0901 

SHIFT 

STORE 

649 

0403 

X 

03 

595 

0715 

PF 

15 

650 

1511 

INDIR 

596 

0601 

RECALL 

651 

0709 

RE 

09 

597 

0800 

SR 

00 

652 

0208 

f 

08 

598 

0609 

ST 

09 

653 

0701 

RE 

01 

599 

0001 

E 1 

6 54 

0607 

ST 

37 

600 

0309 

- 

09 

655 

0001 

E 1 

501 

09  ) 1 

SHIFT 

STORE 

656 

0307 

- 

07 

632 

0930 

SP 

00 

657 

0714 

RE 

14 

0 15 

09  30 

Shift 

MARK 

658 

0207 

07 

604 

0103 

T 

03 

659 

1511 

INDIR 

Figure  14.  Final  Wang  code  listing  (cont'd) . 
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srrp 

COOF 

Kr  V 

STEP 

COHt 

I-JIJT  Tf'N 

Kf"Y 

660 

0707 

PC 

07 

710 

0900 

SHIFT 

MARK 

661 

0607 

ST 

07 

711 

0107 

T 

07 

66;? 

0 70  5 

712 

0100 

T 

CO 

663 

0 40  7 

X 

07 

713 

0702 

PE 

02 

66^ 

0308 

- 

08 

714 

0902 

SHIFT 

ALPHA 

666 

1511 

IND  IP 

715 

C8C6 

S IN 

666 

060  9 

ST 

09 

716 

C80C 

SEARCH 

667 

0702 

Pf 

02 

71  7 

0005 

ALL  UP 

05 

66P 

0600 

S T 

00 

71B 

CE  00 

SEARCH 

669 

0001 

E 1 

719 

00  0 8 

all  up 

08 

670 

0 20  1 

f 

01 

720 

03  00 

- 

00 

6M 

090  2 

SHIFT 

ALPHA 

721 

JlOO 

T 

00 

672 

0806 

SIN 

722 

coo: 

F 0 

673 

0300 

$F APCH 

*72  3 

oooo 

F 0 

674 

0105 

T 

05 

724 

0001 

E 1 

676 

0900 

SHIFT 

MARK 

72  5 

cool 

E 1 

6 76 

0106 

T 

06 

72  6 

0000 

E 0 

677 

0301 

PECALL 

727 

0500 

/ 

on 

679 

0801 

SP 

01 

728 

0200 

4 

00 

679 

0609 

ST 

09 

729 

0100 

T 

00 

68  0 

0711 

pr 

1 1 

730 

0000 

E 0 

681 

0209 

4- 

09 

731 

coco 

B 0 

682 

0 704 

Pb 

04 

732 

CIOC 

T 

00 

683 

151  1 

INDIP 

733 

ClOO 

T 

00 

684 

0609 

ST 

09 

734 

0000 

E 0 

685 

0801 

RECALL 

735 

CIOC 

T 

00 

6 36 

030  1 

SP 

01 

736 

0914 

SHIFT 

END 

68  7 

0609 

ST 

09 

688 

071  2 

FF 

12 

639 

0209 

+ 

09 

690 

0703 

PE 

03 

691 

151  1 

INOTR 

69  2 

0600 

ST 

09 

69  3 

0801 

recall 

694 

0801 

SP 

01 

695 

0609 

ST 

09 

696 

0001 

E I 

697 

0209 

F 

09 

693 

0901 

SHIFT 

STORE 

699 

0801 

SP 

01 

700 

0705 

PF 

05 

701 

0904 

SHIFT 

J NF  0 

70? 

0800 

SEARCH 

703 

0107 

T 

07 

704 

0704 

PE- 

04 

705 

0012 

CHS 

706 

0604 

ST 

04 

707 

0105 

T 

05 

708 

0800 

SEAFCH 

709 

0106 

T 

06 

Figure  14. 

Final  Wang 

code  listing 

(cont * 

d). 
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register  use  must  start  with  a step  number  divisible  by  8,  unused  steps 
preceding  the  data  are  replaced  with  the  GO  code  (0830).  Thus,  there 
could  be  as  many  as  seven  GO  statements  preceding  the  large  number  data 
registers  (the  example  in  fig.  14  has  none) . These  could  be  replaced  by 
a small  user-added  subroutine  or  simply  ignored. 

The  final  printed  output  (fig.  15)  contains  miscellaneous 
information  that  may  be  useful  to  a user  of  the  program.  The  verify 
program  (VP)  number  is  given  for  the  user  to  check  when  the  program  is 
laoded  in  the  machine.  Marks  and  registers  used  are  specified  so  the 
user  can  determine  what  is  left  for  other  purposes.  The  total  number  of 
steps  (program  and  registers)  is  given  for  the  same  purpose.  If  this 
number  is  greater  than  1848,  the  program  is  too  large  for  the  machine. 
Finally,  the  entry  points  are  given  to  facilitate  any  editing  the  user 
may  wish  to  perform. 


ENTRY 

POINTS 

NAME 

MARK 

MAIN 

1000 

POLRT 

1001 

THE  VERIFY  PROGRAM  NUMBER  IS  7832 
REGISTERS  USED 

16  - 133 
MARKS  USED 

0 - 107 

1030  - 1001 

118  REGISTERS  WERE  USED 
1681  STEPS  OF  TOTAL  STORAGE  WERE  USED 
Figure  15.  Entry  point  and  register  information. 
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The  punched  cards  provided  as  output  by  the  compiler  are  put  into  a 
local  program  to  punch  the  actual  Wang  cards  read  by  the  card  reader 
(sect.  2) . The  first  card  shows  the  hexadecimal  number  of  cards 
following,  and  each  succeeding  card  contains  40  steps  at  two  columns  per 
step  in  hexadecimal  representation.  (The  code  0512  is  written  as  5C; 

i.e.,  high-order  and  low-order  code  each  takes  a column.) 

To  use  the  program  once  the  code  is  in  the  machine,  it  is  necessary 
to  know  how  to  control  the  I/O  flow.  Hence,  it  is  wise  to  write  the 
instructions  on  the  back  of  the  final  Wang  input  cards.  For  this 
program,  such  instructions  would  probably  look  like  the  following: 


Purpose:  To  find  roots  of  real  polynomial  ^ ^i+l^^'  0 < M < 25 

i=0  ^ 

(Translation  of  POLRT)  (Explanation) 

VP  = 7719 — load  only  at  step  0000 
Put  Printer  ON. 

1.  Key  fO.  Entry  point  for  MAIN 

2.  4 appears  in  display.  Enter  degree 

M of  polynomial  and  key  GO  (printed  4 is  index  of  M (fig.  6) 
M) 

3.  2.i  appears  in  display.  Enter  C.  2 is  index  of  XCOF 

and  key  GO  (printed  C) . ^ (fig.  6) 

After  last  coefficient  is  entered,  roots  are  computed  and 
printed  (X  real,  I imaginary) , unless  they  are  not  found. 
Then,  error  code  is  printed  (labelled  E) 

E = 1 means  M is  less  than  1.  E=2  should  never  arise 

since  we  have  made 

E = 2 mecins  M is  greater  than  36.  provision  for  M < 26 

E = 3 means  unable  to  determine  root  with  500  iterations  on  5 
starting  values. 

E = 4 means  C , . =0 . 

M+i 

Other  information  could  be  indicated  if  it  might  be  useful  (we  use 
the  symbol  to  mean  "is  stored  in  register")  : 
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M 13 


Coefficients  16  to  40 
Real  part  of  roots  ->^66  to  89 

Corresponding  imaginary  part  of  roots  90  to  113 
Marks  used:  fOO-1,  0000-T07 

Total  steps  used:  1680  (736  for  code  at  low  end  and  944  at  high 

end) 

With  these  instructions,  only  the  source  and  object  listings  need  be 
saved  if  trouble  can  be  anticipated. 

It  is  hoped  that  the  detail  of  this  example  helps  clarify  how  to  use 
the  compiler.  Experiments  will  probably  clarify  the  use  further  and 
suggest  other  techniques  for  improvements.  The  authors  would  appreciate 
comments  or  suggestions  for  this  purpose. 
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APPENDIX  A.— ADAPTABILITY  OF  SYSTEM  TO  OTHER  COMPUTERS 

The  compiler  system  has  been  exclusively  written  in  FORTRAN  IV  and 
is  independent  of  the  word  length  (32  bits)  or  character  code  (EBCDIC) r 
except  for  the  following  cases: 

(a)  FUNCTION  DIG — This  function  takes  a digit  character  code  (left 
justified)  and  right  justifies  as  numerical  digit. 

(b)  FUNCTION  NUMBER“This  function  is  the  opposite  of  DIG. 

The  program  takes  approximately  280k  bytes  of  core  on  the 

IBM  370/195.  In  addition^  two  routines  reference  a direct  access  file 
(logical  unit  1)  that  contains  records  466  words  in  length. 

Approximately  130  records  can  be  stored  on  the  disk.  The  two  routines 
are  listed  below: 


(a) 

GET 

(A, 

I, 

0) 

— Gets 

array 

A 

from 

the 

Ith 

record 

on  the 

file. 

(b) 

PUT 

(A, 

I, 

0) 

— Puts 

array 

A 

into 

the 

Ith 

record 

on  the 

file. 

The  routine  BIN  takes  the  loader  codes  (n  Wang  steps)  and  punches 
out  an  object  deck  containing  40  steps  to  a card^  each  step  represented 
as  a two-digit  hexadecimal  code.  A leader  card  gives  the  count  of  the 
number  of  cards  (also  in  hexadecimal)  in  the  object  deck.  This  object 
deck  is  then  converted  into  the  Wang  card  form  through  a conversion 
routine  written  for  the  IBM  1130.  This  routine  could  be  altered  to 
produce  the  Wang  object  deck  directly. 

A listing  of  the  compiler  is  available  upon  request.  The  system 
contains  approximately  11,000  cards,  including  4000  comment  cards.  The 
source  code  may  be  obtained  by  sending  a tape  to  the  authors.  The  code 
will  be  placed  on  a nine-track  tape  at  800  BPI  on  one  file.  Also 
available  upon  request  is  a listing  of  the  IBM  1130  program  for 
converting  the  object  decks. 

In  converting  the  system  to  a smaller  computer,  it  will  be  necessary 
to  overlay  the  system.  Fortunately,  the  system  is  divided  into  six 
autonomous  subsystems  (sect.  3 of  the  main  body  of  the  report)  with  the 
important  tables  being  passed  through  labelled  commons  or  disk  files. 


65 


DISTRIBUTION 


DEFENSE  DOCUMENTATION  CENTER 
CAMERON  STATION,  BUILDING  5 
ALEXANDRIA,  VA  22314 
ATTN  DDC-TCA  (12  COPIES) 

OFC,  CHIEF  OF  RESEARCH  & DEVELOPMENT 
USA  RSCH  & DEV  GROUP  (EUROPE) 

BOX  15 

FPR  NEW  YORK  09510 
ATTN  LTC  EDWARD  E.  CHICK 

CHIEF,  MATERI7VLS  BRANCH 

COMMANDER 

US  ARMY  MATERIEL  DEVELOPMENT 
& READINESS  COMMAND 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 
ATTN  DRCRD,  RES,  DEV,  & ENGR 
DIRECTORATE 

ATTN  DRCRD-T,  RESEARCH  DIV 
COMMANDER 

USA  ARMAMENT  COMMAND 
ROCK  ISLAND,  IL  61201 
ATTN  DRSAR-ASF,  FUZE  DIV 
ATTN  DRSAR-RDF,  SYS  DEV  DIV  - FUZES 

COMMANDER 

USA  MISSILE  & MUNITIONS  CENTER 
& SCHOOL 

REDSTONE  ARSENAL,  AL  35809 
ATTN  ATSK-CTD-F 

DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
WASHINGTON,  DC  20305 
ATTN  APTL,  DAS A TECH  LIBRARY 

DIRECTOR  OF  DEFENSE  RES  AND 
ENGINEERING 
WASHINGTON,  DC  20301 
ATTN  TECHNICAL  LIBRARY 

DIRECTOR 

NATIONAL  SECURITY  AGENCY 
FORT  GEORGE  G.  MEADE,  MD  20755 
ATTN  T.  A.  PRUGH 

COMMANDER 

US  ARMY  RESEARCH  OFFICE  (DURHAM) 

P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709 
ATTN  CRD-AA-IP 

COMMANDER 

USA  ELECTRONICS  COMMAND 
FORT  MONMOUTH,  NJ  07703 
ATTN  DRSEL-CE,  COMMUNICATIONS- 

ELECTRONICS  INTEGRATION  OFC 
ATTN  DRSEL-TL,  ELECTRONICS  TECHNOLOGY 
& DEVICES  LABORATORY 


USA  ELECTRONICS  COMMAND  (CONT'D) 

ATTN  DRSEL-WL,  ELECTRONIC  WARFARE  LAB 
ATTN  DRSEL-GG,  COMPUTER-AIDED  DESIGN 
& ENGINEERING  OFFICE 
ATTN  DRSEL-GG,  TECHNICAL  LIBRARY 

MOUNTAIN  VIEW  OFFICE  ( DRSEL-WL- RU) 
ELECTRONIC  WARFARE  LABORATORY 
P.O.  BOX  205 
MOUNTAIN  VIEW,  CA  94040 

COMMANDER 

USA  MISSILE  COMMAND 
REDSTONE  ARSENAL,  AL  35809 
ATTN  DRSMI-RBLD,  CHIEF  DOC  SECTION 

COMMANDER 

USA  MOBILITY  EQUIPMENT  R&D  CENTER 
FORT  BELVOIR,  VA  22060 
ATTN  SMEFB-W,  TECHNICAL  LIBRARY 

COMMANDER 
EDGEWOOD  ARSENAL 
EDGEWOOD  ARSENAL,  MD  21010 
ATTN  SMUEA-TS-L,  TECH  LIBRARY 

COMMANDER 
FRANKFORD  ARSENAL 
BRIDGE  & TACONY  STREETS 
PHILADELPHIA,  PA  19137 
ATTN  KIOOO,  TECHNICAL  LIBRARY 

COMMANDER 
PICATINNY  ARSENAL 
DOVER,  NJ  07801 

ATTN  SARPA-TS-T-S , TECHNICAL  LIBRARY 
COMMANDER 

USA  ABERDEEN  PROVING  GROUND 
ABERDEEN  PROVING  GROUND,  MD  21005 
ATTN  STEAP-TL,  TECH  LIBRARY,  BLDG  305 

COMMANDER 

USA  ELECTRONICS  PROVING  GROUND 
FORT  HUACHUCA,  AZ  85613 
ATTN  STEEP-PA-I,  TECH  INFO  CENTER 

COMMANDER 

YUMA  PROVING  GROUND 
YUMA,  AZ  85364 

ATTN  STEYP-MTL,  TEST  ENGINEERING  DIV 
COMMANDER 

USA  WEAPONS  COMMAND,  HA 
ROCK  ISLAND,  IL  61201 
ATTN  SWERR-PL,  TECHNICAL  LIBRARY 

CHIEF  OF  NAVAL  OPERATIONS 
NAVY  DEPARTMENT 
WASHINGTON,  DC  20350 
ATTN  NOP- 09 8,  DIR,  OFC  OF  RES,  DEV, 
TEST,  AND  EVALUATION 
ATTN  NOP-985F,  WEAPONS  TECH  BR 
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NAVAL  ELECTRONICS  LABORATORY  CENTER 
SAN  DIEGO,  CA  92152 
ATTN  TECHNICAL  LIBRARY 

COMMANDER 

PACIFIC  MISSILE  RANGE 
NAVAL  MISSILE  CENTER 
POINT  MUGU,  CA  93042 
ATTN  CODE  5632,  TECHNICAL  LIBRARY 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  MD  20910 
ATTN  L-315,  TECH  LIBRARY 

COMMANDER 

NAVAL  SEA  SYSTEMS  COMMAND 
2521  JEFFERSON  DAVIS  HIGHWAY 
ARLINGTON,  VA  20360 
ATTN  NSEA-0632,  LIBRARY  BRANCH 

DIRECTOR 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  DC  20390 
ATTN  2620,  TECHNICAL  LIBRARY  BR 

COMMANDER 

NAVAL  SHIP  SYSTEMS  COMMAND,  HQ 
2531  JEFFERSON  DAVIS  HIGHWAY 
WASHINGTON,  DC  20360 
ATTN  NSHP-2052,  TECH  LIBRARY  BR 

COMMANDER 

NAVAL  WE7VPONS  CENTER 
CHINA  LAKE,  CA  93555 
ATTN  CODE  753,  LIBRARY  DIV 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN,  VA  22448 
ATTN  TECHNICAL  LIBRARY 

US  AIR  FORCE,  HEADQUARTERS 
DCS,  RESEARCH  & DEVELOPMENT 
WASHINGTON,  DC  20330 

COMMANDER 

HQ  AIR  FORCE  SYSTEMS  COMMAND 
ANDREWS  AFB 
WASHINGTON,  DC  20331 
ATTN  DAPL,  TECHNICAL  LIBRARY 
ATTN  DPSL,  TECH  LIBRARY 

COMMANDER 

AF  CAMBRIDGE  RESEARCH 
LABORATORIES,  AFSC 
L,  G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
ATTN  E.  CZERLINSKY 


COMMANDER 

ARMAMENT  DEVELOPMENT  7VND  TEST  CENTER 
EGLIN  AIR  FORCE  BASE,  FL  32542 
ATTN  ADTC(DLOSL),  TECH  LIBRARY 

COMMANDER 

AERONAUTICAL  SYSTEMS  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
ATTN  ASD/SD,  DEPUTY  FOR  SYSTEMS 
ATTN  TECHNICAL  LIBRARY 

COMMANDER 

HQ  SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
P.  O.  96960  WORLDWAYS  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
ATTN  SN,  DEP  FOR  SPACE  COMM  SYS 
ATTN  SYT,  COMPUTER  TECHNOLOGY  OFC 

COMMANDER 

AF  SPECIAL  WEAPONS  CENTER,  AFSC 
KIRTLAND  AFB,  NM  87117 
ATTN  SWTSX,  SURVIVABILITY/ 

VULNERABILITY  BRANCH 

HQ,  SAAMA,  SANEPA 
KELLEY  AFB,  TX  78241 
ATTN  DIR  OF  MATERIEL  MANAGEMENT 

US  ENERGY  RESEARCH  & DEVELOPMENT 
ADMINISTRATION 
WASHINGTON,  DC  20545 
ATTN  TECHNICAL  LIBRARY 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  DC  20234 
ATTN  LIBRARY 

LIBRARY  OF  CONGRESS 
SCIENCE  & TECHNOLOGY  DIVISION 
WASHINGTON,  DC  20540 
ATTN  HEAD,  LIB  OPNS, 

NASA  AMES  RESEARCH  CENTER 
MOFFETT  FIELD,  CA  94035 
ATTN  S.  J.  DE  FRANCE,  DIRECTOR 

NASA  GEORGE  C.  MARSHALL  SPACE  FLIGHT  CTR 
HUNTSVILLE,  AL  35812 
ATTN  M-G  & C-NS 

NASA  GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
ATTN  LIBRARY 

NASA  LEWIS  RESEARCH  CENTER 
21000  BROOKP7UUC  ROAD 
CLEVELAND,  OH  44135 
ATTN  LIBRARIAN 
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ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13440 
ATTN  LTF,  COMPUTER  ENGINEERING  BR 
ATTN  TECHNICAL  LIBRARY 

NASA  SCIENTIFIC  & TECH  INFO  FACILITY 
P.  0.  BOX  33 
COLLEGE  PARK,  MD  20740 
ATTN  ACQUISITIONS  BR  (S-AK/DL) 

NATIONAL  OCEANIC  & ATMOSPHERIC  ADM 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
BOULDER,  CO  80302 
ATTN  LIBRARY,  R-51,  TECH  REPORTS 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 
JET  PROPULSION  LABORATORY 
4800  OAK  GROVE  DRIVE 
PASADENA,  CA  91103 
ATTN  TDS,  LIBRARY  MANAGER 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  RADIATION  LABORATORY 
BERKLEY,  CA  94720 
ATTN  LIBRARY,  BUILDING  50,  RM  134 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
P.O.  BOX  1663 
LOS  ALAMOS,  NM  87544 
ATTN  R.  GAWLER 

UNIVERSITY  OF  FLORIDA 
GAINSEVILLE,  FL  32603 
ATTN  R.  C.  JOHNSON,  JR. 

ATTN  R.  D.  WALKER 

UNIVERSITY  OF  ILLINOIS 
DEPARTMENT  OF  MATHEMATICS 
URBANA,  IL  61801 
ATTN  LAWRENCE  A WHITE 

UNIVERSITY  OF  MARYLAND 
COMPUTER  SCIENCE  DEPARTMENT 
COLLEGE  PARK,  MD  20741 
ATTN  DR.  YAOHAN  CHU 

UNIVERSITY  OF  MICHIGAN 
INFRARED  INFORMATION 
& ANALYSIS  CENTER 
ANN  ARBOR,  MI  48106 
ATTN  WILLIAM  L.  WOLFE 

BELL  TELEPHONE  LABORATORIES 
WHIPPANY  ROAD 
WHIPPANY,  NJ  07981 
ATTN  LIBRARIAN 


TYMSHARE  INC. 

1911  NORTH  FORT  MYER  DRIVE 
ARLINGTON,  VA  22209 
ATTN  CARLYLE  REEDER 

WANG  LABORATORIES 
8360  NORTH  STREET 
TEWKSBURY,  MA  01876 
ATTN  JASON  TAYLOR 
ATTN  HAROLD  KOPLOW 
ATTN  ROY  KOLK 

ILLINOIS  STATE  WATER  SURVEY 
BOX  232 

URBANA,  IL  61801 
ATTN  MARIE  F.  BURNS,  LIBRARIAN 

DIGITAL  ACOUSTICS,  INC. 

1415  E.  McFADDEN,  SUITE  F 
SANTA  ANA,  CA  92705 
ATTN  MISS  PAMELA  HURST 

COMMANDER 
AFATL/DLRD 
EGLIN  AFB,  FL  32542 
ATTN  MR.  COLLINS 

NATIONAL  INSTITUTES  OF  HEALTH 
BETHESDA,  MD  20014 
ATTN  DR.  C.  PATLAK,  BLDG  13,  RM  1D24 

MGR  SYSTEMS 
1510  RICHARDS  AVENUE 
WILLIAMSPORT,  PA  17701 
ATTN  MR.  BURNETT  TYSON 

CHIEF  ENGR  COM-PU-TOR 
76  OX  YOKE  DRIVE 
WETHERSFIELD,  CT  06109 
ATTN  WM  H . SMYERS , JR. 

THE  JOHNS  HOPKINS  UNIVERSITY 
DEPARIMENT  OF  CHEMISTRY 
BALTIMORE,  MD  21218 
ATTN  DR.  JOYCE  J.  KAUFMAN 
ATTN  MR.  HARRY  J.  T.  PRESTON 

NAVAL  SECURITY  ENGINEERING  FACILITY 
3801  NEBRASKA  AVENUE 
WASHINGTON,  DC  20390 
ATTN  MR.  JOHN  H.  BICKFORD,  CODE  0245 

KINNEY  SHOE  CORPORATION 
16TH  FLOOR 
233  BROADWAY 
NEW  YORK,  NY  10007 
ATTN  MR.  GEORGE  J.  MICHELSON 
SPEC  CONS  TO  THE  PRES 
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ICHTHYOLOGICAL  ASSOCIATES,  INC, 
SCHUYLKILL  RIVER  ECOLOGICAL  STUDY 
FRICKS  LOCK  ROAD,  RD  1 
POTTSTOWN,  PA  19464 
ATTN  MR.  KEN  LITE 

BLECK  ENGINERRING  CO.,  INC. 

1321  GLEN  ROCK  AVENUE 
WAUKEGAN,  IL  60085 
ATTN  DONNA  L.  BLECK 

DIRECTOR  OF  LABORATORIES 
HOLY  CROSS  HOSPITAL 
2701  WEST  68TH  ST 
AT  CALIFORNIA  AVENUE 
CHICAGO,  IL  60629 
ATTN  A.  M.  RING,  MD 

J.  J,  GARCIA  & ASSOCIATES,  INC. 
11039  N,  E,,  6TH  AVENUE 
MIAMI,  FL  33161 
ATTN  EMRIQUE  ALVAREZ 

LINE  COUPLING  EQUIPMENT  ENGINEERING 
GENERAL  ELECTRIC  COMPANY 
MOUNTAIN  VIEW  ROAD 
LYNCHBURG,  VA  24502 
ATTN  MR.  D.  B.  BRAH , MANAGER 

OREGON  STATE  UNIVERSITY 
SCHOOL  OF  OCEANOGRAPHY 
CORVALLIS,  OR  97331 
ATTN  DR.  LOUIS  I.  GORDON 
ASSISTANT  PROFESSOR 

CARTER  PRODUCTS 
RESEARCH  LABORATORY 
CRANBURY,  NJ  08512 
ATTN  W.  M.  WOODING 

DIRECTOR,  TECHNICAL  SERVICES 

PICKARD  & ANDERSON  ENGINEERS 
69  SOUTH  ST 
AUBURN,  NJ  13021 
ATTN  WILLIAM  C.  ANDERSON,  P.E. 

NORTHEASTERN  PRODUCTS  COMPANY 
3500  S.  CLINTON  AVENUE 
SOUTH  PLAINFIELD,  NJ  07080 
ATTN  RICHARD  D.  GUIDO 

MANAGER  QUALITY  CONTROL 

COMMANDANT 

USA  FIELD  ARTILLARY  SCHOOL 
FORT  SILL,  OK  73503 
ATTN  BILL  MILLSPAUGH 
ATSFCTD/SD 


CALCULATOR  CONSULTANT 
45- 3A  MT,  PLEASANT  VILLAGE 
MORRIS  PLAINS,  NJ  07950 
ATTN  MR.  NEAL  H.  KUHN 

US  ENVIRONMENTAL  PROTECTION  AGENCY 
P.O.  BOX  5036 
ROCHESTER,  NY  14627 
ATTN  MR.  DONALD  J.  CASEY 
CHIEF,  IFYGL  BRANCH 

ST.  MARY’S  UNIVERSITY 
2700  CINCINNATI  AVENUE 
SAN  ANTONIO,  TX  28284 
ATTN  DR.  TOM  MOTE 

ETHYL  CORPORATION 
FUNDAMENTAL  STUDIES  DEPARTMENT 
TERRE  HAUTE,  IN  47808 
ATTN  MR.  CHARLES  FURLAND 

BOEING  AEROSPACE  COMPANY 
P.O.  BOX  3999 
SEATTLE,  WA  98124 
ATTN  MR.  MALCOLM  MATHEWS 
MS  8C-41 

PENNWALT  CORPORATION 
TECHNOLOGICAL  CENTER 
900  FIRST  AVENUE 
KING  OF  PRUSSIA,  PA  19406 
ATTN  DR.  J.  E,  DOHANY 

MR.  L.  H.  CHAMBERLIN 
10510  SUNNYBROOK  LANE,  SW 
TACOMA,  WA  98498 

HARRY  DIAMOND  LABORATORIES 
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OFFICER/FLYER,  I. N. /LANDIS,  P.E./ 
SOMMER,  H. /CONRAD,  E.E. 

ATTN  CARTER,  W.W.,  DR.,  ACTING  TECHNICAL 
DIRECTOR/MARCUS,  S.M. 
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ATTN  CHIEF,  DIV  700 
ATTN  CHIEF,  DIV  800 
ATTN  CHIEF,  LAB  900 
ATTN  CHIEF,  LAB  1000 
ATTN  RECORD  COPY,  BR  041 
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